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ABSTRACT 


Spatial  and  temporal  variability  of  the  Black  Sea  surface  circulation  with  the  link  to 
chlorophyll-a  concentration  and  surface  winds  is  investigated  using  Satellite  data  from 
Archiving,  Validation,  and  Interpretation  of  Satellite  Oceanographic  data  (AVISO),  Sea¬ 
viewing  Wide  Field-of-view  Sensor  (SeaWiFS),  and  Quick  Scatterometer  (QuikSCAT) 
with  the  Self-Organizing  Maps  (SOMs).  Six  spatial  patterns  with  temporal  variability 
were  identified  for  the  surface  currents:  Pattern- 1  (Sevastopol  Cyclonic  and  Batumi 
Dipole  Eddies,  21%),  Pattern-2  (Cyclonic  RIM  Current  and  Anti-cyclonic  Batumi  Eddy, 
16%),  Pattern-3  (Anti-cyclonic  Sevastopol  and  Batumi  Eddies,  17%),  Pattern-4  (Cyclonic 
RIM  Current  and  Cyclonic  Batumi  Eddy,  21%),  Pattern-5  (Anti-cyclonic  RIM  Current 
and  Batumi  Dipole  Eddies,  15%),  Pattern-6  (Anti-cyclonic  RIM  Current  and  Multi 
Eddies,  10%).  It  is  found  the  change  of  the  bi-modal  characteristics  in  2000-2009  with 
the  fall  bloom  being  more  significant  than  the  spring  bloom.  The  surface  circulation 
Pattem-4  (cyclonic  RIM  current  and  Batumi  eddy)  is  associated  with  the  occurrence  of 
the  fall  bloom.  Evident  connection  of  negative  NAO  and  negative  ENSO  to  the  Pattern-4 
circulation  implies  the  large-scale  atmospheric  effect.  Possible  connection  of  these 
patterns  to  the  climatological  indices,  such  as  the  North  Atlantic  Oscillation  (NAO)  and 
the  East  Atlantic/West  Russian  (EAWR),  oscillation  are  also  discussed. 
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I.  INTRODUCTION 


A.  THE  BLACK  SEA 

The  Black  Sea  is  located  between  Asia  Minor  and  southeastern  Europe.  It  is 
connected  to  the  Mediterranean  Sea  through  the  Bosporus,  the  Sea  of  Marmara,  the 
Dardanelles  and  the  Aegean  Sea.  If  the  Sea  of  Azov  is  left  out,  the  Black  Sea  covers 
approximately  436,400  km  .  The  Black  Sea  has  been  pursuing  its  crucial  importance  to 
commerce  in  the  region  for  ages  (http://www.ceoe.udel.edu/blacksea/geography/ 
index.html). 

Six  countries  neighbor  the  main  inland  sea.  Namely:  Romania  and  Bulgaria  at  the 
west  coast;  Ukraine  and  Russia  at  the  north  coast;  Georgia  at  the  east  coast;  and  finally 
Turkey  at  the  south  coast.  Apart  from  these  neighboring  countries,  there  are  10  other 
nations  that  affect  the  region  by  the  five  major  rivers,  those  empty  into  the  Black  Sea. 
The  largest  one  among  these  five  rivers  is  the  Danube  River 
(http  ://www  .ceoe  .udel.edu/blacksea/ geography/index  .html) . 


Figure  1.  The  Black  Sea  and  the  surrounding  countries 
(Map  available  from  Geography. Org,  http://world-geography.org/sea/92-black-sea.html) 
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The  strategic  importance  of  the  Black  Sea  region  has  been  increasing 
significantly,  especially  since  the  memberships  of  Bulgaria  and  Romania  were  granted 
for  North  Atlantic  Treaty  Organization  (NATO)  in  2004,  and  shortly  after  for  the 
European  Union  (EU)  in  2007.  Due  to  the  fact  that  The  Black  Sea  became  the  border 
region  of  both  NATO  and  the  EU,  any  issues  affecting  the  geo-economic,  geo-political  or 
the  geo-strategic  status  of  the  region  draw  attention  of  both  establishments  effectively. 
Examples  of  the  current  geo-political  and  geo-strategic  issues  are:  the  frozen  conflicts  of 
Moldova  (Transnistria)  and  Georgia  (Abkhazia-South  Ossetia),  both  of  which  are 
fundamentally  based  on  the  ethnicity  diversity;  Armenia- Azerbaijan  (Nagorno-Karabakh) 
which  arose  as  a  result  of  the  tortious  occupancy  of  Azerbaijan  soil  by  the  Armenian 
forces;  the  resurge  of  Russia;  organized  criminal  activity  and  transnational  terrorism  (J. 
Bosbotinis  2013,  unpublished  manuscript).  Therefore,  it  is  clear  that  expansion  of  the 
borders  of  both  NATO  and  EU  confronts  both  communities  with  the  problems  and 
weaknesses  of  the  Black  Sea  area  as  well  as  its  benefits,  and  challenges  the  Euro-Atlantic 
goals  in  terms  of  keeping  the  continent  secure,  free  and  wealthy  (Ulger  2007). 


Legend 


|Moldova 
j  Transnistria 
I  lother  states 


Figure  2.  The  map  showing  the  current  situation  in  Moldova  (Transnistria)  conflict 

(From  Wikimedia  Commons ) 
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The  geo-economic  importance  of  the  Black  Sea  is  expected  to  rise  as  oil  and  gas 
pipelines  from  Azerbaijan  and  Russia  centers  in  the  region.  Its  importance  may  be  even 
more  pronounced,  if  these  resources  are  integrated  with  a  Pan-European  transport  system. 
The  increasing  trend  of  the  gas  demand  is  expected  to  continue  in  the  following  years. 
Especially,  the  estimate  of  the  gas  consumption  in  Europe  is  at  least  600  billion  cubic 
meters  after  2020,  while  domestic  production  will  drop  down  to  half.  Therefore,  South 
Stream  is  determined  as  a  feasible  project  that  will  help  European  countries  to  overcome 
the  crisis  of  the  lack  of  energy  sources  in  the  horizon,  and  at  the  same  time  it  will 
conserve  the  energy  security  (http://www.gazprom.com/press/news 
/2010/november/article  106074/).  Furthermore,  the  Black  Sea  region  neighbors  Caspian 
Basin,  which  is  an  extremely  important  region  due  to  its  hydrocarbon  resources.  This 
proximity  increases  strategic  importance  of  the  Black  Sea  as  the  energy  security  becomes 
the  focus  of  both  NATO  and  the  EU  (J.  Bosbotinis  2013,  unpublished  manuscript). 
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South  Stream  pipeline  route  options 

-  Route  option  1 :  Russia  -  Bulgaria  -  Serbia  -  Hungary  -  Austria 

—  Route  option  2:  Russia  -  Bulgaria  -  Serbia  -  Hungary  -  Slovenia 
Route  option  3:  Russia  -  Bulgaria  -  Serbia  -  Hungary  -  Austria 
and  Bulgaria  -  Greece  -  Italy 
““  Route  common  to  all  the  above  options 


Figure  3.  An  illustration  showing  the  options  for  the  South  Stream  pipeline  route  (From 

Gazprom) 
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Russia  is  the  other  actor  in  the  same  sense,  and  tries  to  maintain  its  dominance  in  the 
wider  Eurasian  post-Soviet  space  and  monopolistic  status  in  terms  of  energy  exports. 
Russia  also  pays  extreme  attention  to  the  region  and  its  issues.  Namely:  the  futures  of  the 
North-South  Caucasus  and  Ukraine;  expanding  borders  of  NATO  and  the  EU;  the 
increasing  trend  of  US  presence.  As  a  result  of  the  common  interests  in  the  Black  Sea,  the 
region  evolves  into  a  core  area  where  a  competition  between  U.S.  and  Russia  takes  place. 
Both  actors  try  to  take  advantage  of  the  Black  Sea  and  its  surrounding  by  treating  it  as  a 
gateway  to  extend  their  influence;  particularly  Middle  East,  Central  Asia  and  South 
Caucasus  areas  are  of  interest  for  the  US;  whereas,  Mediterranean  in  addition  to  these 
regions  for  Russia  (J.  Bosbotinis  2013,  unpublished  manuscript). 

B.  OBJECTIVE 

The  new  militarily,  politic  and  economic  developments  in  the  Black  Sea  propels 
the  need  of  a  deeper  understanding  the  nature  and  the  mechanisms  of  the  basin.  Hence, 
the  spatial  and  temporal  variation  of  the  Black  Sea  holds  a  crucial  importance,  especially 
with  regard  to  operational  purposes.  The  main  focus  of  this  research  is  to  obtain  data  sets 
concerning  the  three  features  of  the  Black  Sea  (geostrophic  velocity  anomalies, 
geostrophic  winds,  and  chlorophyll-a  concentration)  through  various  systems,  and  to 
analyze  them  to  extract  this  spatial  and  temporal  variation.  To  identify  a  correlation 
among  these  patterns  of  separate  data  sets  and  to  evidence  that  one  pattern  of  a  data  set 
can  affect  another  from  a  different  data  set  is  intended  afterwards.  There  are  number  of 
tools  available  to  achieve  this  goal.  The  self-organizing  map  (SOM)  is  selected 
particularly  in  this  study,  due  to  the  fact  that  its  accuracy  is  greatly  appreciated  and  it  is 
not  widely  used  in  the  Black  Sea  region. 

C.  SELF-ORGANIZING  MAP 

Even  though  data  availability  advances  steadily,  both  by  in-situ  observations  and 
remote-sensing  techniques,  it  doesn’t  necessarily  correspond  to  the  efficient  and 
consistent  processing  or  usage  of  these  data.  The  estimation  is  that  no  more  than  5%  of 
the  images  collected  by  remote- sensing  techniques  are  either  exposed  to  human  eyes  or 
actually  processed.  Reasons  for  this  may  include  lack  of  time  or,  from  a  different 
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perspective,  lack  of  tools  that  can  process  a  huge  amount  of  data  efficiently,  while 
extracting  the  key  features  that  this  thesis  is  after  (Petrou  2004). 

The  SOM  is  introduced  as  an  option  to  abridge  this  ever  widening  gap  between 
the  available  and  processed  data.  SOM  (aka,  the  Kohonen  map)  is  an  unsupervised  neural 
network  and  competitive  learning  underlies  it  (Kohonen  1998,  2001).  One  of  the  key 
features  of  SOM  is  its  capability  of  preserving  the  neighborhood  relations  of  the  high¬ 
dimensional  input  data  while  projecting  it  onto  a  lower  dimensional  (preferably  two- 
dimensional)  space.  Therefore,  it  is  known  as  a  topology-preserving  technique.  As  a 
result  of  this  topology-preserving  technique,  patterns  that  are  similar  are  put  into 
neighboring  regions  during  the  mapping  process;  having  these  similar  patterns  closely 
mapped  to  one  another  enables  a  huge  advantage  in  pattern  recognition  and  the  feature 
extraction  processes  of  the  complex  data  sets.  The  SOM  and  its  benefits  were  presented 
to  the  oceanography  community  by  Richardson  et  al.  (2003),  since  then  it  has  been 
heavily  used  in  miscellaneous  projects.  Some  of  these  works  are  listed  in  Table  1  (Liu 
and  Weisberg  2011). 
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Table  1.  A  few  examples  of  the  SOM  applications  in  oceanography  (From  Liu  and 

Weisberg  2011). 


Oceanographic  data 

Regions 

References 

Oil  spill 

Galician  coast, 

Europe 

Corchado  et  al.  (2008),  Mata  et  al. 
(2009), 

Borges  et  al.  (2010) 

Maritime  data 

Europe 

Lobo  (2009) 

Salinity 

Columbia  River 

plume 

Liu  et  al.  (2009) 

Surface  winds 

Southeast  Atlantic 

Richardson  et  al.  (2003),  Risien  et  al. 
(2004) 

Satellite  measured 

sea  surface  height 

Southeast  Atlantic 

Hardman-Mountford  et  al.  (2003) 

Indian  Ocean 

Iskandar  (2009) 

South  China  Sea 

Liu  et  al.  (2008) 

Satellite  ocean  color, 
chlorophyll 

Pacific 

Ainsworth  (1999),  Ainsworth  and  Jones 
(1999) 

Southeast  Atlantic 

Yacoub  et  al.  (2001) 

Southwest  Atlantic 

Saraceno  et  al.  (2006) 

North  Atlantic 

Niang  et  al.  (2003),  Telszewski  et  al. 
(2009) 

D.  THESIS  ORGANIZATION 

Chapter  II  provides  a  foundation  for  the  Black  Sea  dynamics.  Atmospheric 
forcing  mechanism,  surface  currents,  temperature-salinity  properties,  and  biological 
environment  are  touched  on. 

Chapter  III  gives  information  about  the  data  used  here  in  this  research,  where  the 
data  was  obtained,  how  the  data  was  collected,  as  well  as  a  look  at  the  technology  behind 
the  data  collection. 

Chapter  IV  details  the  SOM  in  terms  of  foundation.  The  idea  behind  this 
technique,  how  it  works,  the  mathematics  that  it  is  based  on,  and  the  toolbox  for  Matlab 
will  be  mentioned  briefly. 
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Chapter  V  presents  the  analyses  that  were  applied  to  the  data  set  of  the 
geostrophic  velocity,  as  well  as  the  results  of  these  analyses.  The  relation  between  the 
geostrophic  velocity  and  the  large-scale  phenomena  are  discussed. 

Chapter  VI  presents  the  analyses  and  their  results  regarding  the  geostrophic  wind 
data  set.  As  in  Chapter  V,  the  relation  between  the  geostrophic  velocity  and  the  large- 
scale  phenomena  are  discussed. 

Chapter  VII  focuses  on  the  chlorophyll-a  data  set  and  investigates  the  correlation 
between  chlorophyll-a  and  geostrophic  velocity,  chlorophyll-a  and  geostrophic  wind 
separately. 

Chapter  VIII,  the  final  chapter,  concludes  the  research  by  pointing  out  the 
findings  and  manifests  any  recommendations  for  further  study. 
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II.  THE  DYNAMICS  OF  THE  BLACK  SEA 


A.  INTRODUCTION 

The  Black  Sea  is  a  mid-latitude  basin  located  between  41°N  and  46°N.  It 
exchanges  water  with  the  Mediterranean  Sea  through  the  Strait  of  Bosphorus.  In  terms  of 
bathymetry,  the  Black  Sea  has  a  flat  abyssal  plain  at  the  center,  which  comprises  the 
maximum  depth  found  (2200  m);  and  along  its  coast  there  is  a  continental  shelf  that  has 
various  offshore  extensions  depending  on  the  location.  The  extension  off  the  western 
Turkish  coast,  for  instance,  is  ~5  km,  whereas  it  is  -200  km  off  the  northwestern  coast. 
The  Black  Sea  is  split  into  two  sub-basins  by  the  Crimean  Peninsula  to  the  north  and  the 
central  Turkish  coastline  to  the  south.  Apart  from  the  Mediterranean  Sea,  the  Black  Sea  is 
also  fed  by  large  river  discharges,  specifically  by  the  Danube,  the  Dnepr,  and  the  Dnestr 
Rivers  which  are  located  on  the  northwestern  shelf.  Therefore,  the  Black  Sea  houses  both 
the  high  salinity  waters  of  Mediterranean  origin  and  the  low  salinity  waters  of  riverine 
origin.  The  high  salinity  waters  underlay  the  low  salinity  waters.  This  formation  results  in 
a  permanent  pycnocline,  which  is  formed  around  100-150  m  beneath  the  surface  and 
prevents  any  exchanges  between  these  two  waters  of  different  salinity  profiles.  Deeper 
waters  that  are  below  500  m  are  immobile.  During  the  summer  time  and  between  10-40 
m,  formation  of  a  seasonal  thermocline  is  observed  (Oguz  et  al.  1992). 

B.  SURFACE  CURRENTS 

The  dominant  characteristic  of  the  Black  Sea  circulation  is  the  basin-wide 
circulation  that  is  mainly  cyclonic.  This  circulation  depends  heavily  on  time  and  is 
considered  as  following  the  continental  slope.  This  unique  characteristic  is  called  the  Rim 
Current  and  it  is  40-80  km  wide.  It  is  fundamentally  fueled  by  the  mean  cyclonic  wind 
pattern  in  the  area  and  by  the  intense  buoyancy  input  (Stanev  1990;  Oguz  et  al.  1995; 
Korotaev  et  al.  2001).  According  to  the  model  simulations,  apart  from  the  forcing 
mechanism  of  the  winds,  the  Rim  Current  disappears  unless  the  topology  data  is  put  into 
the  model  (Oguz  et  al.  1995).  Furthermore,  if  the  bottom  slope  and  the  course  of  the 
coastline  are  changed  along  the  Turkish  and  Caucasian  coasts,  the  Rim  Current  meanders 
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and  features  of  instability  are  widely  observed,  in  terms  of  both  time  and  space  (Oguz  et 
al.  1993).  Models  have  also  demonstrated  that  the  weakening  of  the  wind  stress  in  the 
spring  and  summer  results  in  a  drop  in  the  mean  current’s  intensity  (Gregoire  et  al.  2004; 
Stanev  1990).  This  drop  in  the  intensity  causes  intensified  meandering  of  the  Rim  Current 
along  the  Turkish  and  Caucasian  coasts  resulting  in  the  creation  of  large  (100-200  km) 
meanders.  The  presence  of  a  series  of  recurrent,  near-shore,  anti-cyclonic  eddies  between 
the  Rim  Current  and  the  coast,  along  with  a  number  of  cyclonic  gyres  in  the  basin’s 
central  area,  have  been  confirmed  by  both  satellite  data  (Oguz  et  al.  1992;  Sur  et  al.  1994; 
Stanev  et  al.  2000;  Ginzburg  et  al.  2000;  Korotaev  et  al.  2001,  2003)  and  by  hydrographic 
observations  (Oguz  et  al.  1993,  1994;  Poulain  et  al.  2005). 


Among  these  near-shore  and  anti-cyclonic  eddies,  two  persistent  eddies  are 
noticeable;  namely:  the  Batumi  Eddy  and  the  Sevastopol  Eddy.  The  Batumi  Eddy  is 
located  off  of  the  southeastern  coast  and  the  Sevastopol  Eddy  is  located  west  of  the 
Crimean  Peninsula  and  over  the  gentle  continental  slope.  The  interior  part  of  the  Rim 
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Current  is  composed  of  one  of  two  detached  cyclonic  gyres.  These  detached  cyclonic 
gyres  consist  of  either  one  gyre  located  in  the  eastern  part  of  the  sea  and  the  other  located 
in  the  western  part  of  the  sea,  or  a  single  extended  basin-wide  cell.  At  the  times  it  is 
composed  of  two  separate  cyclonic  gyres,  recurrent  mesoscale  eddies  form  in  the  interior 
parts  of  the  gyres  and  they  are  connected  to  one  another  via  a  recurrent  anticyclone 
known  as  the  Central  Basin  Eddy  (Oguz  et  al.  1993).  The  scales  of  the  central  gyres  and 
eddies  vary  from  tens  to  hundreds  of  kilometers.  Additionally,  the  number,  location  and 
morphology  of  the  central  gyres  and  eddies  are  strongly  season-dependent  (Poulain  et  al. 
2005).  Eddies  that  vary  inter-annually,  seasonally,  and  by  mesoscale/synoptic  scale, 
dominate  the  near-surface  circulation  of  the  Black  Sea  (Oguz  et  al.  1994;  Poulain  et  al. 
2005). 

According  to  Zatsepin  et  al.  (2003),  apart  from  the  Rim  Current  zones  and  the 
near-shore,  anti-cyclonic  eddies  also  form  in  the  central  Black  Sea.  There  is  an 
interaction  between  these  eddies  and  the  Rim  Current,  which  causes  mixing  of  near-shore 
waters  and  advection  due  to  reflection  from  the  coast  and  the  two  separate  jets  that  it 
forms.  Additionally,  there  is  no  correlation  between  the  anticyclones  in  the  eastern  part  of 
the  central  Black  Sea  that  have  long  life  periods  and  the  harshness  of  the  preceding 
winter.  On  the  other  hand,  the  variability  of  the  wind  is  a  decisive  factor  affecting  the 
scales  of  synoptic,  seasonal,  and  inter-annual  eddies.  Exchanges  in  the  shelf-deep  basin 
of  the  Black  Sea  are  mainly  driven  by  the  eddy  dynamics  (Zatsepin  et  al.  2003;  Poulain  et 
al.  2005). 

The  Black  Sea  eddies  are  formed  in  the  eastern  part  of  the  basin  and  they 
propagate  to  the  western  part  of  the  basin  in  the  form  of  Rossby  Waves  at  the  speed  of 
~3cm/s.  This  eddy  propagation  is  strongly  affected  by  the  narrow  part  of  the  Black  Sea 
located  at  the  south  of  the  Crimean  Peninsula  (Rachev  and  Stanev  1997).  The  more 
eddies  lose  their  momentum  and  shrink  to  smaller  scales  in  the  western  basin,  the  more  a 
rise  in  the  dissipation  is  observed.  This  process  is  influenced  by  topography.  Korotaev  et 
al.  (2001)  examined  sea  level  data  that  was  collected  by  satellite  altimeters  and  confirmed 
this  westward  propagation  of  phase  (Poulain  et  al.  2005). 
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Inter-annual  variability  of  the  Black  Sea  circulation  is  significantly  affected  by 
river  discharges,  local  dynamics,  and  wind  forcing’s  seasonal  variability  (Stanev  et  al. 
1995;  Stanev  and  Beckers  1999).  Mesoscale,  seasonal,  and  inter-annual  variability  was 
best  described  by  the  observations  made  by  satellite  altimeters  (Stanev  et  al.  2000; 
Korotaev  et  al.  2001,  2003).  According  to  Stanev  et  al.  (2000),  intensification  of  the 
general  cyclonic  circulations  occurs  in  spring  and  winter,  depending  on  the  significant 
increase  in  sea  level  in  the  near-shore  areas.  The  wind  forcing  is  also  pointed  out  as  the 
major  factor  that  drives  the  circulation’s  intensification  at  the  seasonal  scale.  The  Rim 
Current  and  large  coastal  quasi-permanent,  anti-cyclonic  eddies  (the  Batumi  Eddy,  the 
Sevastopol  Eddy,  etc.)  have  the  strongest  variations  at  the  intra-annual  scale.  The  Batumi 
Eddy  shows  clear  seasonal  variability  in  terms  of  circulation  such  that  it  demonstrates  a 
pronounced  cyclonic  circulation  during  the  winter  and  a  weaker  cyclonic  or  an  anti- 
cyclonic  circulation  during  the  summer  and  fall.  Unlike  the  Batumi  Eddy  case,  the 
variability  of  the  Sevastopol  Eddy  does  not  depend  on  seasons,  but  rather  follows  an 
inter-annual  or  intra-annual  scale  (Poulain  et  al.  2005). 

Korotaev  et  al.  (2001)  used  both  TOPEX-POSEIDON  and  ERS-1  data  to  explain 
that,  there  is  a  strong  variability  of  the  Black  Sea  circulation  due  to  seasonality.  More 
specifically,  the  circulation  intensifies  during  the  winter  and  spring,  and  it  weakens 
during  the  summer  and  fall.  The  period  of  the  most  intense  mesoscale  oscillations  is 
determined  to  be  120  days.  These  oscillations  are  found  in  the  basin’s  southeastern  part. 
The  Rim  Current  splits  and  meanders  off  the  Crimean  Peninsula  in  this  part  of  the  basin 
due  to  the  influence  of  the  local  bathymetry  (Poulain  et  al.  2005). 

The  general  circulation  of  the  Black  Sea  is  comprised  of  two-cyclonic  gyres 
during  the  winter,  which  convert  into  a  single  cell  structure  which  is  encompassed  by  a 
weaker  and  wider  Rim  Current  during  the  summer.  The  Rim  Current  demonstrates  an 
intense  and  complicated  mesoscale  activity  during  the  summer  and  fall  when  the  whole 
circulation  structure  decomposes  into  a  number  of  interconnecting  eddies  (Korotaev  et  al. 
2003;  Poulain  et  al.  2005).  The  dominance  of  the  anti-cyclonic  gyres  in  the  northwestern 
and  southeastern  corners  of  the  Black  Sea  during  the  summer  and  fall  is  confirmed  by  the 
model  simulations  as  well.  Additionally,  the  simulations  noted  that  the  most  pronounced 
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and  persistent  coastal  anti-cyclonic  eddy  is  the  Batumi  Eddy,  which  mostly  dominates 
throughout  the  period  between  March  and  October.  These  simulations  also  showed  that 
the  conditions  during  the  winter  and  summer  months  favor  the  formation  of  the 
Sevastopol  Eddy  (Poulain  et  al.  2005). 

C.  WIND  FIELD 

The  meteorological  observations  concerning  the  Black  Sea  region  give  clues 
about  the  high  variability  and  distinction  of  the  circulation  patterns  of  the  atmospheric 
processes.  The  main  external  factor  that  induces  the  drift  currents  is  considered  as  wind 
stress.  Accompanied  by  coastline,  baroclinicity,  and  bottom  relief,  wind  stress  constructs 
a  particular  Black  Sea  circulation  under  the  influence  of  Earth’s  rotation  (Kordzadze  et  al. 
2007). 


Figure  5.  Illustration  of  the  temporal  variability  of  the  wind  stress  from  the  5 -minute- 
resolution  Black  Sea  Modular  Ocean  Model  (MOM)  (Stanev  et  al.  2003),  along  with  the 
matching  curve  calculated  from  the  European  Centre  for  Medium-Range  Weather 
Forecasts  (ECMWF) -reanalyzed  data  (From  Stanev  2005). 

Black  Sea  oceanography  aims  to  investigate  the  evolution  of  Black  Sea 
circulation  while  taking  the  atmospheric  circulation  strongly  into  account.  For  this 
purpose,  the  quantitative  circulation  aspects  of  the  Black  Sea  should  be  put  into  use  as 
well  as  the  experimental  methods  (Kordzadze  et  al.  2007). 
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The  Black  Sea’s  east  coast  is  the  region  where  some  notable  local  wind  fields  are 
located  due  to  the  coastal  mountain  range  bordering  the  sea  that  has  irregular  height  with 
gaps  and  valleys.  Easterly  air  flows  to  the  Black  Sea  through  these  gaps  and  valleys.  In 
the  Black  Sea’s  southern  region,  a  broad  valley  called  Kolkhida  enables  this  airflow  to 
progress  even  more  westward.  The  Novorossiyskaya  bora,  which  can  exceed  the  speed  of 
40  m/s,  is  the  most  well-known  local  east  coast  wind  of  the  Black  Sea  (Gusev  1959; 
Burman  1969;  Alpers  et  al.  2009).  This  local  wind  can  be  very  dangerous,  particularly  for 
coastal  traffic  and  those  participating  in  harbor  activities  (Kordzadze  et  al.  2007). 

1.  BORA  WINDS 

When  a  cold  front  system  passes  over  a  region  or  when  a  high  pressure  gradient 
occurs  and  moves  the  cold  air  over  a  coastal  mountain  range,  local  down-slope  winds  are 
observed  and  these  winds  are  called  Bora  winds  (Alpers  et  al.  2010).  These  winds  are 
found  in  mountainous  coastal  regions  with  heights  below  600  m,  so  that  descending  cold 
air  warms  adiabatically  at  a  relatively  small  rate.  This  is  the  direct  opposite  condition 
when  compared  to  Foehn  and  Chinook  winds,  where  the  adiabatic  warming  rate  is 
notable.  These  winds  are  also  known  as  Novorossiyskaya  Boras  due  to  the  fact  that  they 
are  mostly  observed  in  the  coastal  region  close  to  Novorossiysk  (Gusev  1959;  Burman 
1969).  At  times  these  winds  with  notable  strength  stay  over  the  region  for  several  days, 
they  excite  a  cyclonic  eddy  at  the  southeast  of  the  Tuapse  in  Russia  (Alpers  et  al.  2010). 
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Figure  6.  A  satellite  image  demonstrating  a  Bora  event  taking  place  over  the  east  coast 
of  the  Black  Sea  on  December  15,  2008  at  19:10  UTC  (From  ESA) 

2.  KATABATIC  WINDS 

Katabatic  winds  blow  following  a  sloping  coastal  topography  and  eventually  they 
reach  the  surface  of  the  sea.  Besides  being  cold  winds.  Katabatic  winds  are  observed  only 
in  the  early  morning,  in  the  evening,  or  at  night  (Alpers  et  al.  2010).  The  inducing 
mechanism  for  these  types  of  winds  is  the  cooling  rate  difference  between  air  near  the 
surface  over  the  land  and  over  the  sea.  Since  air  near  the  surface  has  a  larger  cooling  rate 
over  the  land  as  opposed  to  over  the  sea,  the  cold  air  flows  downhill.  This  flow  is  known 
as  gravity  flow.  Unless  low-level  clouds  inhibit  the  cooling  process,  air  is  able  to  cool  off 
notably  fast  during  the  evening  and  night.  This  fast  cooling  of  the  air  enables  the 
generation  of  the  Katabatic  winds.  Since  the  wind  flowing  downhill  is  channeled  through 
coastal  valleys,  the  katabatic  winds  form  a  tongue-like  shape  along  the  coastline  (Alpers 
et  al.  1998;  Alpers  et  al.  2010). 
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Figure  7.  A  satellite  image  demonstrating  the  sea  surface  signatures  of  a  Katabatic  wind 
event  taking  place  over  the  east  coast  of  the  Black  Sea  on  April  30,  2008  at  19:07  UTC 

(From  ESA) 


3.  FOEHN  WINDS 

The  driving  mechanism  for  the  warm  Foehn  winds  is  descending  motion  of  the 
airflow.  Therefore,  they  are  observed  on  the  lee  side  of  mountain  ranges.  They  are 
accompanied  by  ascending  temperatures,  diffusion  of  low-level  clouds,  and  descending 
relative  humidity.  Apart  from  the  eastern  part  of  the  Black  Sea,  Foehn  winds  are  mostly 
observed  on  the  north  side  of  the  Alps  (Alpers  et  al.  2010). 
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Figure  8.  A  satellite  image  demonstrating  the  sea  surface  signatures  of  a  Foehn  wind 
event  taking  place  over  the  southeast  coast  of  the  Black  Sea  on  January  11,  2010  at  07:31 

UTC  (From  ESA) 
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Figure  9.  Monthly  averaged  climatological  wind  fields.  The  arrow  scale  equals  1  m/s 

(From  Peneva  and  Stips  2004) 


D.  TEMPERATURE  AND  SALINITY 


As  a  result  of  the  abundance  of  freshwater  sources  (e.g.,  rivers  and  rainfalls) 
compared  to  lost  freshwater  due  to  evaporation,  the  Black  Sea  holds  a  positive  freshwater 
balance.  More  specifically,  the  evaporation  removes  around  354  square  km  of  fresh 
water,  whereas  freshwater  inputs  from  river  outflows  and  precipitations  are 
approximately  350  square  km  and  occur  approximately  230  times  each  year,  respectively. 
This  positive  freshwater  balance  of  the  Black  Sea  results  in  an  average  of  0.43m  gain  in 
sea  level  in  comparison  with  the  Marmara  Sea. 

The  excessive  amount  of  water  in  favor  of  the  Black  Sea  runs  into  the  Marmara 
Sea  via  the  strait  of  Bosphorus.  Therefore  the  straits  comprise  two  different  flows.  One  of 
them  is  the  upper  flow,  which  transports  the  surface  water  from  the  Black  Sea,  and  the 
other  is  the  bottom  flow  that  transports  35%o  of  the  Mediterranean  origin  salt  water  to  the 
Black  Sea.  The  salinity  of  the  surface  waters  is  17.5  -  18%o.  It  is  relatively  low  due  to  the 
mixing  process  of  the  salty  waters  with  the  original  waters  of  the  basin.  In  terms  of  the 
volume,  the  bottom  flow  is  approximately  300  cubic  km,  which  is  around  half  of  the 
surface  flow. 

The  surface  salinity  has  seasonal  variability.  In  terms  of  the  vertical  distribution 
of  salinity,  an  increase  is  observed  at  50  m  depth;  when  the  depth  hits  200  m  the  increase 
rate  is  still  remarkable,  whereas  below  200  m  the  increase  rate  is  much  slower.  The  salt 
composition  of  the  waters  of  the  Black  Sea  resembles  the  composition  of  oceans.  The 
Black  Sea  waters  house  remarkable  measures  of  nutrients,  especially  the  mixtures  of 
nitrogen  and  phosphorus  that  originate  in  the  rivers.  Even  though  the  surface  waters’ 
temperature  shows  variability  throughout  the  year,  this  is  not  the  case  for  the  water 
column  below  the  surface.  For  these  deeper  waters,  a  rather  unique  vertical  change  in 
temperature  is  observed.  The  depths  that  the  minimum  temperatures  have  been  observed 
start  from  50-60  m  go  down  to  80-90  m.  For  the  deeper  parts  of  the  water  column,  a 
slow  increase  in  the  temperature  is  noted  until  reaching  the  seabed  (~  2200  m). 

The  density  of  the  seawater  is  assessed  by  the  vertical  distribution  of  salinity  and 
temperature  values.  Due  to  the  unique  features  the  Black  Sea  introduces,  it  is  comprised 
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of  two  diverse  water  layers:  the  first  measuring  from  0  to  about  200  m  (called  the  lighter 
upper  layer),  and  the  second  layer  measuring  between  200  m  and  the  seabed  (  called  the 
heavier  lower  layer).  The  weak  vertical  circulation  arises  within  the  waters  due  to  this 
stratification.  The  mixing  of  these  distinct  layers  does  not  occur  efficiently.  Furthermore, 
life  in  the  Black  Sea  is  significantly  influenced  by  this  structure 
(http://www.blackseascene.net/content/content. asp?menu=0040032_000000). 
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Figure  10.  Monthly  averaged  mixed-layer  temperature  of  the  Black  Sea  obtained  through  the  climatological  run  of  General  Estuarine 

Transport  Model  (GETM)  (From  Peneva  and  Stips  2004) 
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Monthly  averaged  mixed-layer  salinity  of  the  Black  Sea  obtained  through  the  climatological  run  of  General  Estuarine 

Transport  Model  (GETM)  (From  Peneva  and  Stips  2004) 
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Figure  12.  Illustration  of  the  temporal  variability  of  the  thermal  buoyancy  flux  from  the 
5-minute-resolution  Black  Sea  Modular  Ocean  Model  (MOM)  (Stanev  et  al.  2003),  along 
with  the  matching  curve  calculated  from  the  European  Centre  for  Medium-Range 
Weather  Forecasts  (ECMWF)-reanalyzed  data 
(From  Stanev  2005). 
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Figure  13.  Illustration  of  the  temporal  variability  of  the  haline  buoyancy  flux  from  the  5- 
minute -resolution  Black  Sea  Modular  Ocean  Model  (MOM)  (From  Stanev  2005). 


According  to  model  simulations  (Stanev  and  Staneva  2001),  the  peak  of  the 
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water’s  heating  (2  x  10  W/m")  occurs  in  June,  whereas  the  peak  of  the  cooling  (3  x  10 

W/m")  takes  place  in  early  winter.  Moreover,  the  amplitude  of  the  haline  buoyancy  flux  is 
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surpassed  by  the  amplitude  of  the  thermal  buoyancy  flux  by  an  order  of  magnitude.  On 
the  other  hand,  in  terms  of  net  values,  the  net  haline  buoyancy  flux  is  approximately  four 
times  larger  than  the  net  thermal  buoyancy  flux. 

E.  CHLOROPHYLL-A 

The  concentration  of  chlorophyll-a  is  the  most  suitable  and  useful  parameter  in 
terms  of  monitoring  the  marine  ecosystem.  For  instance,  there  is  a  high  correlation 
between  the  level  of  the  chlorophyll-a  concentration  and  the  primary  production;  that 
level  is  also  determined  as  the  indicator  of  the  phytoplankton  biomass.  Since  the 
concentration  of  chlorophyll-a  has  been  monitored  for  many  years  and  still  continues, 
there  are  many  data  sets  available  for  research  purposes  (Yunev  et  al.  1987;  Berseneva 
1993;  Vedernikov  and  Demidov  1993;  Yilmaz  et  al  1998;  Yunev  et  al  2002;  Kopelevich 
et  al.  2002). 

Satellite  ocean-color  data  provides  the  synoptic  scale  measurements  of  marine 
ecosystems,  and  the  concentration  of  chlorophyll-a  is  an  exclusive  feature  of  these 
ecosystems.  Temporal  and  spatial  distributions  of  the  chlorophyll-a  concentration  are 
identified  by  satellite  ocean-color  sensors  lately,  which  are  widely  accepted  as  the 
standard  tool  for  this  particular  data  collection  (Kopelevich  et  al.  2002).  The  Coastal 
Zone  Color  Scanner  (CZCS)  was  designed  as  the  first-generation  ocean-color  sensor  and 
Nimbus-7,  a  NASA  satellite,  was  the  carrier  of  this  sensor.  During  the  period  of  October 
1978  to  June  1986,  a  large  amount  of  beneficial  information  regarding  ocean  color  was 
collected  via  the  CZCS.  This  period  is  considered  the  stable  phase  of  the  Black  Sea 
ecosystem.  The  severe  changes  in  the  Black  Sea  ecosystem  began  shortly  after  that 
period.  These  changes  were  accompanied  by  the  rush  of  the  “Mnemiopsis  leidyi,”  which 
was  carried  to  the  region  mistakenly  since  the  summer  of  1988  (Shuskina  and 
Vinogradov  1991).  These  changes  make  the  CZCS  measurements  even  more  important 
since  they  allowed  researchers  to  identify  the  period  of  stable  conditions  in  the  Black  Sea, 
and  these  changes  enabled  researchers  to  compare  current  conditions  with  the  stable 
conditions  in  order  to  better  understand  the  subsequent  occurrences  (Kopelevich  et  al. 
2002).  Due  to  the  CZCS’  lack  of  significant  coverage  of  the  Black  Sea  region  during  cold 
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periods,  it  is  not  feasible  to  study  the  interannual  fluctuations  between  warm  and  cold 
periods  in  the  region  (Nezlin  et  al.  1999;  Kopelevich  et  al.  2002). 

According  to  the  mean  monthly  distributions  provided  by  Kopelevich  et  al. 
(2002),  the  distinction  between  the  regions  of  the  western  shelf  and  the  Black  Sea’s  open 
region  is  obvious.  The  eastern  shelf,  on  the  other  hand,  appears  to  be  confined  adjacent  to 
the  coastline,  while  the  regions  of  the  western  shelf  demonstrate  great  variability  in  terms 
of  their  mean  monthly  distributions.  The  open  regions  show  relatively  uniform  patterns, 
especially  for  the  period  between  April  and  October  (Kopelevich  et  al.  2002). 
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Figure  14.  The  average  monthly  distributions  of  chloropyll-a  (mg/m  )  found  in  the  CZCS  data  for  the  first 

6-month  period  (From  Kopelevich  et  al.  2002). 
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Figure  15.  The  average  monthly  distributions  of  chlorophyll-a  (mg/m  )  found  in  the  CZCS  data  for  the  second 

6-month  period  (From  Kopelevich  et  al.  2002). 


27 


In  1993,  Vedernikov  and  Demidov  studied  the  surface  chlorophyll-a 
concentrations  based  on  different  parts  of  the  basin  during  different  seasons.  They  chose 
November  through  March  as  the  cold  period,  May  through  September  as  the  warm 
period,  and  the  months  of  April  and  October  as  the  transition  periods.  According  to  this 
study,  the  mean  chlorophyll-a  concentrations  at  the  surface  of  the  open  parts  of  the  Black 
Sea  were  found  to  be  0.97-1.52  mg/m3  during  the  cold  period,  0.28-0.38  mg/m3  during  the 
warm  period,  0.61  mg/m  during  April,  and  finally  0.43  mg/m  during  October.  These 
results  are  well-correlated  with  the  outcomes  of  the  study  done  by  Kopelevich  et  al.  2002 
which  selected  the  same  periods.  In  the  latter  study,  the  concentrations  found  for  the  open 
parts  of  the  basin  were  appeared  to  be  0.70-0.83  mg/m  during  the  cold  period,  0.14-0.15 

Q  TO 

mg/m  during  the  warm  period,  0.19-0.20  mg/m  during  April,  and  0.18-0.25  mg/m 
during  October.  Both  studies  detected  the  lowest  concentrations  during  warm  periods  and 
the  highest  concentrations  during  the  cold  periods.  The  differences  between  the 
concentrations  of  chlorophyll-a  are  justified  by  the  additional  time  frame  in  the  latter 
study  (Kopelevich  et  al.  2002). 

The  open  parts  of  the  western  shelf  and  the  eastern  shelf  seemed  to  share  the  same 
mean  chlorophyll-a  concentration  in  terms  of  mean  annual  and  warm  period  values. 
There  was  an  obvious  discrepancy  for  the  cold  period  as  the  mean  concentration  for  the 
open  parts  of  the  eastern  shelf  was  found  to  be  approximately  19%  more  in  comparison 
with  that  of  the  western  shelf.  During  the  transition  months,  however,  the  western  shelf 
demonstrated  a  10-28%  more  concentration  (Kopelevich  et  al.  2002). 

The  coastal  parts  of  the  eastern  shelf  exhibited  higher  mean  concentration  values 
compared  to  the  open  parts  throughout  the  year.  This  difference  was  around  40%  during 
the  cold  period,  more  than  twice  the  amount  during  the  warm  period,  and  finally  80% 
higher  for  the  annual  mean.  In  terms  of  the  seasonal  changes,  on  the  other  hand,  both  the 
coastal  and  the  open  parts  of  the  Black  Sea  had  the  same  nature.  During  the  cold  period, 
the  chlorophyll-a  concentrations  were  detected  as  four  to  six  times  more  than  the  amount 
during  the  warm  period  and  three  to  five  times  more  than  the  amount  during  April  and 
October  (Kopelevich  et  al.  2002). 
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For  the  western  shelf,  the  mean  chlorophyll-a  concentrations  found  in  the  open 
parts  and  the  coastal  parts  of  the  Black  Sea  differed  siginificantly,  particularly  for  the 
warm  period.  The  Romanian  coast,  which  is  heavily  influenced  by  the  Danube  River 
flow,  was  found  to  have  the  highest  mean  chlorophyll-a  concentrations  throughout  the 
entire  year.  During  the  warm  period,  this  region’s  waters  had  values  18  times  more  as 
compared  with  the  open  parts  and  approximately  nine  times  more  as  compared  to  the 
eastern  shelf.  Since  the  homogeneity  in  the  water  column  rises  during  the  cold  period,  the 
discrepancies  among  the  different  regions  (west  shelf  vs.  east  shelf,  open  region  vs. 
coastal  region,  etc.)  of  the  basin  become  less  pronounced.  Even  though  the  whole  Black 
Sea  basin  is  found  to  be  the  most  homogeneous  during  the  cold  period,  open  parts, 
particularly,  tend  to  exhibit  the  highest  homogeneity  during  the  warm  period  (Kopelevich 
et  al.  2002). 

In  previous  studies,  the  assessment  of  the  long  term  variability  of  the  Black  Sea 
physical  properties  is  made  using  either  satellite  data  (e.g.  Stanev  and  Peneva,  2002; 
Ginzburg  et  al.,  2004;  Nardelli  et  al.,  2010;  Grayek  et  al.,  2010)  or  in-situ  data  (e.g.  Oguz 
et  al.,  2006;  Oguz  and  Gilbert,  2007).  In  addition  to  that,  impacts  of  the  North  Atlantic 
Oscillation  (NAO),  the  East  Atlantic/West  Russian  (EAWR)  (e.g.  Oguz  et  al.,  2006)  and 
the  El  Nino  Southern  Oscillation  (ENSO)  on  the  Black  Sea  hydrodynamics  were 
investigated  (e.g.  Gingzburg  et  al.,  2010;  Oguz  et  al.,  2006).  Aside  from  those  studies, 
Arthur  et  al.  (2012)  quantify  the  correlations  between  the  variability  of  different  physical 
processes  concerning  the  Black  Sea  internal  dynamics  in  their  study,  using  the  SOM  for 
the  wind  and  EOF  for  the  Sea  Level  Anomaly  (SLA)  with  connection  to  Sea  Surface 
Temperature  (SST). 

The  main  aim  of  this  research  can  be  summarized  as  follows:  (1)  to  utilize  an 
SOM  to  detect  recurrent  patterns  of  the  surface  winds  for  the  10-year  period  that  follows 
the  study  of  Arthur  et  al.  (2012),  (2)  to  utilize  an  SOM  to  detect  recurent  patterns  of  the 
surface  currents  for  the  same  10-year  period,  (3)  to  investigate  the  influence  of  the  large- 
scale  teleconnection  indices  on  recurrent  patterns  of  both  data  sets,  (4)  to  correlate 
recurrent  patterns  of  the  both  data  sets  with  each  other,  (5)  to  superimpose  the 
chlorophyll-a  concetration  on  recurrent  patterns  of  both  data  sets. 


29 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


30 


III.  DATA  SOURCES 


A.  AVISO 

Geostrophic  velocity  anomalies  used  in  this  research  were  obtained  through  the 
Archiving,  Validation,  and  Interpretation  of  Satellite  Oceanographic  data  (AVISO), 
which  is  a  system  located  at  Collecte  Localisation  Satellites  (CLS).  CLS  was  founded  in 
1996  and  the  specializations  of  this  processing  center  are  data  collection,  trajectory  of 
satellites,  as  well  as  satellite -based  ocean  observations.  Conductivity-Temperature- 
Depth-Oxygen  Profiler  (CTDP)  is  another  system  run  by  both  CLS  and 
American  Physical  Oceanography-Distributed  Active  Archive  Center  (Po-Daac).  AVISO 
is  connected  to  CTDP,  which  is  capable  of  downloading  Geophysical  Data  Records 
(GDR),  Interim  Geophysical  Data  Records  (IGDR)  and  processed  data  from  both  centers. 
On  the  other  hand,  French  geophysical  corrections,  orbit  files  for  the  French,  and 
American  instruments  are  downloaded  by  AVISO  (http://www.aviso.oceanobs.com 
/no_cache/en/newsstand/newsletter/  newsletterOl/what-is-aviso/). 

1.  ALTIMETRY 

a.  History 

The  discussion  of  using  radar  instrumentation  in  space  for  oceanography 
purposes  first  took  place  in  1969  at  the  congress  of  Williamstown.  Thus,  measurement  of 
the  sea  level  using  radar  techniques  was  initially  aimed  by  satellite  altimetry.  First 
altimetry  measurements  tried  to  determine  the  topography  of  the  ocean  surface.  With  the 
missions  Skylab  and  Geos-3,  the  United  States  became  the  first  country  that  ever  used  the 
satellite-borne  altimeter  technology  in  history.  Following  the  United  States,  in  1991,  the 
European  Space  Agency  (ESA)  launched  ERS-1  and  got  involved  to  this  concept. 

Centre  National  d’Etudes  Spatiales  (CNES)  was  working  on  placing  an 
altimeter  called  Poseidon  onboard  of  Spot-1  in  1981.  Concurrently,  NASA  was 
evaluating  the  Topex  (Topography  Experiment).  Besides  sharing  interest  in  new  missions 
they  had  one  more  detail  in  common:  both  agencies  were  suffering  from  small  budgets 

that  would  not  allow  them  to  execute  the  missions  by  themselves.  In  1983,  negotiations 
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began  regarding  how  to  combine  these  two  projects  together,  and  in  1987  the  partnership 
was  officially  established.  Instruments  were  going  to  be  launched  on  a  United  States 
satellite  using  a  European  rocket,  and  the  project  was  named  Topex/Poseidon  (T/P). 

Topex/Poseidon  was  critical  for  the  oceanography  community,  as  it 
provided  increased  knowledge  of  oceanography  and  improved  the  community’s 
understanding  in  terms  of  ocean  observation.  Even  the  in-situ  observations  made  over  the 
previous  one  hundred  years  were  not  capable  of  providing  the  same  amount  of 
information  that  was  collected  during  the  ten-day  cycle  of  Topex/Poseidon  measurements 
(http://www.aviso.oceanobs.com/en/altimetry.html). 


Figure  16.  A  Topex/Poseidon  cycle  compared  to  100  years  of  in  situ  data  (From 
http://www.aviso.oceanobs.com/en/altimetry/history.html) 


32 


b.  Principle 

High-frequency  signals  (at  least  1700  pulses  Hz)  are  emitted  to  Earth  by 
the  altimeters  and,  once  they  are  reflected,  the  echo  from  the  sea  surface  is  received  by 
the  altimeters  again.  The  basic  purpose  of  this  exercise  is  to  find  out  the  satellite-surface 
distance.  The  time  it  takes  for  a  radar  pulse  to  travel  from  the  satellite  to  the  surface  and 
return  to  the  satellite  again  after  reflection  is  measured.  Additionally,  it  is  known  that 
electromagnetic  waves  travel  at  the  speed  of  light.  So,  the  distance  is  calculated  in  light 
of  this  information.  The  presence  of  the  atmosphere  complicates  this  process  since  the 
water  vapor  that  it  contains,  or  the  ionization  taking  place  in  it  can  decelerate  the 
electromagnetic  waves.  Thus,  the  results  should  be  adjusted  according  to  the 
interferences  of  the  atmosphere.  There  is  potential  for  extracting  other  kinds  of 
information  through  altimetry.  One  of  them  being  the  characteristics  of  the  surface  that 
reflection  of  the  signal  takes  place  and  it  can  be  determined  by  focusing  on  the  shape  and 
the  magnitude  of  the  echoes. 

Tracking  altimetry  satellites  is  another  crucial  process  and  should  be  done 
accurately.  There  are  a  variety  of  ways  for  tracking  satellites,  one  of  which  is  a  system 
called  Doris.  Doris  is  a  location  system  that  was  developed  by  ONES,  and  put  onboard  of 
Topex/Poseidon.  It  is  connected  to  a  network  of  50  transmitting  ground  beacons  spread 
worldwide.  The  Doppler  shift  on  the  beacon  signals  is  used  by  Doris  to  identify  the 
satellite’s  velocity  as  it  orbits  (http://www.aviso.oceanobs.com/en/altimetry.html). 
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Figure  17.  The  concepts  and  systems  that  contribute  to  the  altimetry  measurement  (From 

Cnes) 


c.  Multi-Satellites 

There  should  be  a  reasonable  balance  between  the  temporal  and  spatial 
resolutions  of  an  altimetry  satellite  when  determining  its  orbit.  One  way  to  overcome  this 
compromise  is  to  use  Multi-Satellites  with  different  orbits/resolutions,  and  to  take 
advantage  of  the  multiple  views  of  the  same  spot.  This  process  enables  fine  tuning  of  the 
separate  measurements  (http://www.aviso.oceanobs.com/en/altimetry.html). 


Figure  18.  Comparison  of  the  details  of  images  collected  from  a  single  satellite  and 

multi-satellites  (From  CLS/Cnes) 
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2. 


MISSIONS 


Since  it  was  first  founded,  the  system  AVISO  has  been  conducting  satellite 
altimetry  missions  to  provide  an  immense  international  user  community  with  data  that  is 
collected  through  precise  and  high-quality  altimetry  measurements  with  the  help  of 
advanced  technology  and  better  orbit  determination.  For  instance,  Topex/Poseidon  (T/P) 
had  been  the  pioneer  in  highly  precise  altimetry  measurements  and  it  changed  the 
perspective  significantly  in  terms  of  systematical  sea  surface  height  observations. 

There  are  ongoing  missions,  while  some  have  been  terminated.  Future  missions 
are  planned  to  start  in  2013.  Certain  missions  are  executed  by  standalone  agencies, 
whereas  others  have  the  contribution  of  groups.  The  U.S.  Navy,  the  European  Space 
Agency  (ESA),  the  National  Aeronautics  and  Space  Administration  (NASA),  Centre 
National  d’Etudes  Spatiales  (CNES),  the  China  Academy  of  Space  Technology,  the 
Indian  Space  Research  Organization  (ISRO),  the  European  Organisation  for  the 
Exploitation  of  METeorological  SATellites  (EUMETSAT),  and  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  are  agencies  that  have  been  involved  in  at  least  one 
of  the  missions  conducted. 


Table  2.  Specifications  of  the  Satellites  Used  in  the  Past  Missions  of  AVISO 


Satellite 

Geosat 

ERS-1 

T/P 

GFO 

ERS-2 

Envisat 

Launch 

Date 

03/10/1985 

07/17/1991 

08/10/1992 

02/10/1998 

04/21/1995 

03/01/2002 

End  Date 

01/31/1990 

03/31/2000 

01/18/2006 

1 1/26/2008 

07/06/2011 

May  2012 

Altitude 

800  km 

785  km 

1336  km 

800  km 

785  km 

782.4- 
799.8  km 

Inclination 

N/A 

98.52° 

66° 

O 

oo 

o 

98.52  0 

98.55  0 

Repetitivity 

N/A 

33,35,  168 
days 

9.9156  days 

17  days 

35  days 

30-35  days 

Agency 

U.S.  Navy 

ESA 

NASA  and 
CNES 

U.S.  Navy 
and  NOAA 

ESA 

ESA 

Goal 

Describe 
the  marine 
geoid 

Observe 

Earth  and  its 
environment 

Measure  sea 

surface 

height 

Measure 

ocean 

topography 

Observe 

Earth  and  its 
environment 

Observe 
Earth’s 
atmosphere 
and  surface 
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GEOS  3  SEASAT  GEOSAT  ERS-1  T/P  JASON  ENVISAT  ALTIKA 


Figure  19.  Altimetry  measurement  accuracy  evolution  since  the  first  mission  (From 

Cnes) 


Table  3.  Specifications  of  the  Satellites  being  Used  in  the  Current  Missions  of 

AVISO 


Satellite 

HY-2 

Cryosat 

Jason- 1 

Jason-2 

Launch  Date 

08/15/2011 

04/08/2010 

12/07/2001 

06/20/2008 

Altitude 

971  km 

717  km 

1336  km 

1336  km 

Inclination 

99.3° 

92° 

on 

ON 

o 

ON 

ON 

o 

Repetitivity 

14,  168 
days 

369  days 

9.9156  days 

9.9156  days 

Agency 

China 
Academy 
of  Space 
Technology 

ESA 

CNES  and 
NASA 

CNES, 

NASA, 
EUMETSAT 
and  NO  A  A 

Goals 

Observe 
the  ocean 
dynamics 

Polar 

observation 

Measure  sea 

surface 

height 

Measure  sea 

surface 

height 
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The  Spot  satellite  family  began  their  mission  with  the  very  first  member  Spot-1, 
which  was  in  service  from  February  22,  1986  to  2001.  The  last  two  members  of  the 
family  are  still  active  and  working.  Spot-4  was  launched  on  March  24,  1998,  and  Spot-5 
was  launched  on  May  4,  2002.  The  altitude  for  Spot  family  satellites  is  set  at  830  km.  The 
agency  behind  them  is  CNES  and  their  goal  is  to  observe  Earth. 


Table  4.  Specifications  of  the  Satellites  Planned  to  be  Used  in  the  Future  Missions  of 

AVISO 


Satellite 

Saral 

Sentinel-3 

Jason-3 

Jason-CS 

Swot 

Launch  Date 

2013 

2014 

April  2014 

End  of  2017 

2020 

Altitude 

-800  km 

814.5  km 

1336  km 

1336  km 

970  km 

Inclination 

98.55  0 

98. .65° 

66° 

66° 

78° 

Repetitivity 

35  days 

27  days 

10  days 

10  days 

22  days 

Agency 

ISRO  and 
CNES 

ESA 

CNES, 

NASA, 
EUMETSAT 
and  NO  A  A 

CNES, 

NASA, 
EUMETSAT 
and  NO  A  A 

CNES 

and 

NASA 

Goals 

Observe 
the  oceans 

Deliver 

routine 

operational 

services 

Measure  sea 

surface 

height 

Measure  sea 

surface 

height 

Cover 

more 

than 

90%  of 

the 

globe 

Future  missions  aim  for  better  temporal  and  spatial  resolutions  to  enable  studies 
that  phenomena  and  variations  are  more  closely  investigated.  In  particular,  Saral  will  be 
equipped  with  Ka-band  which  is  going  to  provide  better  observations  of  ice,  rain,  wave 
heights,  and  coastal  zones.  On  the  other  hand,  the  altimeter  on  Sentinel-3  will  feature 
both  low  and  high  resolution  modes.  The  low  resolution  mode  will  be  used  over  oceans 
and  other  homogenous  surfaces,  whereas  the  high  resolution  mode  will  be  used  over  sea 
ice  and  for  other  applications  that  require  more  detail 
(http://www.aviso.oceanobs.com/en/missions.html). 
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3. 


DATA 


There  are  five  different  products  available:  namely  sea  surface  height,  surface 
geostrophic  velocity  (or  anomaly),  wind  and  wave,  auxiliary  data,  indicators,  and  in-situ 
products.  There  are  various  options  for  each  of  these  products,  as  well.  The  users  can 
work  on  either  a  global  or  regional  scale;  they  can  use  real  time  (RT),  near  real  time 
(NRT),  or  delayed  time  (DT)  versions  of  these  products;  and  users  even  can  reach  the 
data  from  multi-satellites  which  is  called  “merged  data.” 

The  delayed  time  product  comes  with  two  versions:  Reference  (Ref)  and  Updated 
(Upd).  Reference  series  are  produced  by  two  satellites  at  a  given  time  (Jason  2-Envisat, 
Jason  1-Envisat,  or  Topex/Poseidon-ERS)  and  they  are  homogenous.  Whereas  Updated 
series  are  fed  by  up  to  four  satellites  ( Jason- 2-Jason-  1-Envisat  since  2009,  or  Jason- 1- 
Topex/Poseidon-Envisat-GFO  from  October  2002  to  September  2005),  and  they  are 
treated  by  a  better  verification  process,  yet  they  are  inhomogeneous  in  terms  of  quality. 

Products  are  distributed  by  the  partnership  of  Segment  Sol  multimissions 
d'ALTimetrie,  d'Orbitographie  et  de  localisation  precise  (Ssalto),  and  the  Data 
Unification  and  Altimeter  Combination  System  (Duacs).  Both  systems  process  the  data 
that  comes  from  all  the  altimeter  missions  and  they  aim  to  provide  the  users  with 
consistent  information  that  is  as  homogenous  as  possible.  Ssalto,  additionally,  operates 
facilities  on  earth  in  order  to  control  the  Poseidon  and  Doris  instruments. 

The  product  used  in  this  research  is  delayed  time  (Upd)  and  merged  mean  sea 
level  anomalies  (DT-MSLA)  for  the  Black  Sea  region.  The  content  of  this  product  is  sea 
surface  height  on  a  regular  grid  with  the  resolution  of  1/8°  x  1/8°,  calculated  with  respect 
to  a  mean  sea  level  profile  as  well  as  geostrophic  velocity  anomalies  that  correspond  to 
height.  The  data  is  available  weekly,  and  provided  each  Wednesday  of  the  week  that  the 
data  is  averaged  over  (http://www.aviso.oceanobs.com/en/data/products/sea-surface- 
height-products/regional/msla-black-sea.html;  “The  altimeter  products  were  produced  by 
Ssalto/Duacs  and  distributed  by  Aviso,  with  support  from  Cnes 
(http  ://w  w  w .  aviso  .oceanobs  .com/duacs/)”) . 
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B.  SEAWIFS 


The  Sea-viewing  Wide  Field-of-View  Sensor  (SeaWiFS)  Project  aims  to  collect, 
archive,  and  distribute  validated  research  data  regarding  ocean  bio-optical  properties  on  a 
global  scale  and  gathered  by  an  ocean  color  sensor  orbiting  around  the  Earth.  Different 
quantities  and  types  of  microscopic  marine  plants  (marine  phytoplankton)  affect  the 
ocean  color  in  a  variety  of  ways.  More  specifically,  the  higher  the  concentration  of 
chlorophyll  and  plant  pigments  other  than  chlorophyll,  the  greener  the  water  color.  The 
goal  here  is  the  derivation  of  that  concentration,  using  ocean  observations  made  by  a 
satellite  and  evaluating  the  ocean  color  from  those  observations 
(http://oceancolor.gsfc.nasa.gov/seawifs/background/seawifs_background.html). 
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Figure  20.  An  example  of  the  images  collected  by  SeaWiFS  Project  (From 

NASA/GSFC) 
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1. 


OBJECTIVES 


The  goal  of  SeaWiFS  data  is  to  study  the  influence  of  oceanic  factors  on  global 
change.  Understanding  the  role  of  oceans  in  the  biogeochemical  cycles  (e.g.,  carbon  cycle 
etc.),  which  is  studied  by  using  the  chlorophyll’s  and  other  marine  phytoplankton 
productions’  magnitude  and  variability,  is  intended.  For  that  purpose,  the  scientific 
community  has  been  using  SeaWiFS  data  to  identify  the  limits  of  highly  abundant  growth 
(spring  blooms),  as  well  as  the  length  of  time  over  which  this  growth  occurs. 


Figure  21.  A  photograph  of  the  scanner  assembly  used  in  SeaWiFS  (From 
http://oceancolor.gsfc.nasa.gov/seawifs/background/figures/brochure_figl.gif) 


The  objective  of  the  project  has  been  to  accurately  collect  ocean  color  data  on  a 
global  scale  for  at  least  5  years  and  extract  reasonable  biological  parameters 
(photosynthesis  rates,  etc.)  from  this  data  with  the  help  of  auxiliary  information.  The  next 
step  has  been  to  make  these  findings  available  to  researchers  around  the  world 
(http://oceancolor.gsfc.nasa.gov/seawifs/background/seawifs_970_brochure.html). 
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2. 


ORGANIZATION 


SeaWiFS  data  is  broadcasted  by  OrbView-2  (OV-2)  satellite  in  real  time  to  the 
Goddard  Space  Flight  Center  (GSFC),  the  High-Resolution  Picture  Transmission  (HRPT) 
station,  and  other  stations,  whereas  Global  Area  Coverage  (GAC)  and  Local  Area 
Coverage  (LAC),  recorded  data  in  other  words,  is  broadcast  to  the  GSFC,  the 
ORBIMAGE,  and  the  WallopsFlight  Facility  (WFF).  Next,  HRPT,  GAC,  LAC,  and 
HRPT  data  are  delivered  to  the  SeaWiFS  Data  Processing  System  (SDPS).  The  SDPS  is 
comprised  of  GSFC  mission  operations;  SeaWiFS  data  processing;  data  capture;  and 
calibration,  validation,  and  data  quality  elements. 


Major  Elements  of  At- Sea  Data  Collection  |  Major  Elements  of  SDPS 


Figure  22.  Major  elements  of  SeaWiFS  Project  (From 
http://oceancolor.gsfc.nasa.gov/seawifs/background/figures/brochure_fig3.gif) 

Raw  data  from  the  satellite  is  received  by  the  data  processing  element  and 

standard  global  ocean  color  data  products  are  generated  out  of  that  raw  data.  The  GSFC 
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Earth  Observing  System  Data  and  Information  System  (EOSDIS)  Distributed  Active 
Archive  Center  (DAAC)  is  responsible  for  archiving  SeaWiFS  data  and  distributing  it  to 
researchers.  Then,  it  receives  that  standard  global  ocean  color  data. 

The  calibration  and  validation  element  establishes  the  procedures  of  calibration 
for  SeaWiFS  data  and  updates  those  procedures.  In  order  to  do  that,  this  element  collects 
auxiliary  data  (i.e.,  wind,  atmospheric  pressure,  etc.),  analyzes  trends  and  variations,  and 
validates  higher-level  products  with  correlative  data. 

The  mission  planning  element  links  SeaWiFS  Project  Office  (SPO)  and 
ORBIMAGE  together  to  enable  a  connection  that  tracks  automatic  data  measurements 
and  transmissions,  schedules  the  scientific  maintenances,  and  catches  problems. 
Additionally,  the  orbit  elements  and  navigation  information  are  provided  by  this  element 
to  some  extent.  Mission  management,  the  space  segment,  and  the  command  and  data 
acquisition  station  are  ORBIMAGE’s  responsibilities  (http://oceancolor.gsfc.nasa.gov 
/seawifs/background/seawifs_970_brochure). 

3.  DATA  COLLECTION 

Two  different  scientific  data  are  produced  by  SeaWiFS  in  terms  of  spatial 
resolutions:  LAC  and  GAC.  Even  though  the  broadcast  for  LAC  is  continuous,  it  is 
selectively  recorded.  The  recording  process  for  GAC,  on  the  other  hand,  is  continuous 
and  is  executed  on  board  of  the  satellite.  Direct  broadcasts  of  LAC  for  the  East  coast  of 
the  United  States  have  been  received  by  the  GSFC  and  other  real-time  broadcasts  of  LAC 
have  been  received  by  a  number  of  global  HRPT  stations.  Recorded  data  of  GAC  and 
LAC  are  first  received  by  NASA  WFF,  and  then  forwarded  to  SDPS  at  GSFC. 

Acquiring  full  GAC  data  every  48  hours  is  essential  for  the  science  goals  of  the 
project.  Therefore,  almost  all  of  the  recording  capacity  onboard  of  the  spacecraft  and 
transmission  time  to  GSFC  is  assigned  for  that  requirement.  Limited  space  left  from  the 
GAC  is  used  for  LAC  recording.  The  LAC  recording  process  has  been  allotted  in  priority 
order:  1)  watching  critical  sensor  functions,  2)  conducting  essential  optical  calibration 
and  validation  processes,  and  3)  collecting  data  for  science  studies  that  are  in  need  of  data 
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in  full  resolution  (http://oceancolor.gsfc.nasa.gov/seawifs/background 

/seawif  s_97  O_brochure  .html) . 

4.  SPACECRAFT  AND  SENSOR 

SeaWiFS  instruments  are  carried  by  the  SeaStar  spacecraft  which  was  built  by 
Orbital  Sciences  Corporation  and  was  launched  on  August  1,  1997  to  Earth’s  low  orbit. 
The  control  and  command  of  the  satellite  are  undertaken  by  the  company  GeoEye.  After 
20  days  following  the  launch,  the  satellite  reached  its  final  orbit  which  is  705  km  circular, 
noon,  sun-synchronous.  Data  collection  was  enabled  30  days  after  the  launch.  Orbit 
determination  is  backed  up  by  a  redundant  amount  of  global  positioning  system  (GPS) 
receivers. 

The  SeaWiFS  instrument  is  comprised  of  two  parts:  an  optical  scanner  and  an 
electronic  module.  The  scanning  mechanisms  operate  a  spinning  half-angle  mirror  and  an 
off-axis  folded  telescope.  The  telescope  spins  at  twice  the  speed  of  the  mirror  and  the 
mirror  is  phase-synchronized  with  the  telescope. 

The  telescope,  accompanied  with  the  half-angle  scan  mirror  setting,  enables  a  design 
configuration  that  keeps  the  level  of  polarization  at  minimum  beyond  58.3  degrees  which 
is  the  maximum  scan  angle  necessity. 

The  signals  those  are  detected,  and  amplified,  are  directed  to  the  electronic 
module  from  the  scanner.  This  electronic  module  is  the  element  where  the  signals  are 
amplified  even  more  strongly  and  then  filtered  for  the  purpose  of  noise  bandwidth 
limitation  (http://oceancolor.gsfc.nasa.gov/seawifs/seastar/spacecraft.html). 
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L-101 «  Ai'cratl  Launch 


Third  Stage  Coast 


Second  Stage  Burnout 


Third  Stago  Ignition 


Third  Stage  Burnout 
and  Orbital  Insertion 


Figure  23.  Launch  diagram  of  SeaWiFS  (From 
http://oceancolor.gsfc.nasa.gov/seawifs/seastar/seawifs_launch_diagram.jpeg) 


Table  5.  Mission  Characteristics  of  SeaWiFS 


Orbit  Type 

Sun  Synchronous  at  705  km 

Equator  Crossing 

Noon+20  min,  descending 

Orbital  Period 

99  min. 

Swath  Width  - 1 

2,801  km  LAC/HRPT  (58.3  °) 

Swath  Width  -  2 

1,502  km  GAC  (45  °) 

Spatial  Resolution 

1.1  km  LAC,  4.5  km  GAC 

Revisit  Time 

1  day 
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C.  QUIKSCAT 

1.  SCATTEROMETRY 
a.  Introduction 

Scatterometers  are  airborne  or  space-borne  microwave  radar  sensors  that 
are  designed  to  transmit  a  pulse  of  a  microwave  signal  to  the  surface  of  the  earth  in  order 
to  measure  the  echoed  microwave  energy  while  scanning. 

The  main  research  area  that  space-bome  scatterometry  is  used  for  has  been 
primarily  the  near-surface  winds  over  the  ocean.  These  instruments  are  specifically  called 
the  wind  scatterometers.  They  provide  radar  cross-sections  of  the  surface  accompanied  by 
views  of  different  azimuth  angles,  which  help  researchers  better  understand  the  near¬ 
surface  wind  vector  that  is  present  over  the  ocean  surface  when  processed  by  a 
geophysical  model  function  (GMF).  GMFs  provide  information  regarding  the 
relationship  between  wind  and  backscattered  energy.  The  radar  echoes  coming  from  the 
ocean  surface  are  produced  by  wind- generated  capillary-gravity  waves  which  mostly 
correspond  to  the  near-surface  wind  over  the  ocean.  This  scattering  mechanism  is  called 
Bragg  scattering  (named  after  physicist  Sir  William  Lawrence  Bragg),  and  it  results  from 
waves  that  correspond  with  the  microwaves. 


Figure  24.  Principle  of  Scatterometer  (From  NOAA) 
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Wind  speed  and  direction  determines  the  backscattered  power.  This 
backscattered  power  from  waves  varies  when  observed  from  different  azimuth  angles. 
These  variations  lead  to  the  process  of  the  estimation  and  determination  of  sea  surface 
wind  speeds,  which  is  known  as  wind  retrieval. 

Data  obtained  using  ocean  scatterometers  are  crucial  to  the  research 
community,  especially  in  the  area  of  ocean  circulation  and  air-sea  interaction,  and  their 
influence  on  climate  and  weather  patterns.  Ocean  scatterometers  also  enable  the  research 
of  unusual  phenomena,  such  as  La  Nina  or  changes  that  take  place  in  polar  regions  over 
the  long  term.  These  data  have  also  been  used  for  the  verification  and  justification  of  the 
outputs  that  are  produced  by  the  atmospherical,  oceanic  and  coupled  models,  to  improve 
modeling  accuracy  and  increase  the  forecast  period 
(http://winds.jpl.nasa.gov/aboutScat/index.cfm). 

b.  History 

Information  regarding  the  winds  above  the  ocean  surface  used  to  come 
from  irregular  and  mostly  inaccurate  reports  from  ships  and  buoys. 

Scatterometry  originated  from  early  radar  technology  used  in  World  War 
II.  The  measurements  made  by  this  early  radar  were  found  to  be  noisy,  but  that  noise  was 
discovered  to  be  the  signal  response  echoed  from  the  surface  winds.  This  relationship 
between  the  radar  response  and  surface  winds  wasn’t  established  until  late  1960s. 

Skylab  missions  in  1973  and  1974  were  the  first  examples  of  space-borne 
scatterometry  and  they  convinced  the  community  that  it  is,  in  fact,  beneficial  and  feasible 
to  conduct  scatterometry.  Following  those  missions,  the  Seasat-A  Satellite  Scatterometer 
(SASS)  was  put  into  service  in  October  1978  and  demonstrated  that  wind  velocity  can  be 
measured  from  space  accurately.  European  Space  Agency's  Remote  Sensing  Satellite 
(ERS-1)  was  the  next  mission  in  that  category  and  a  single-swath  scatterometer  was 
onboard  of  the  spacecraft. 

The  NASA  Scatterometer  (NSCAT)  had  the  privilege  of  housing  the  first 
dual-swath  Ku-band.  It  was  launched  in  August  1996  onboard  of  a  Japanese  Advanced 
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Earth  Observing  Satellite  (ADEOS).  Between  the  activation  of  the  instrument  in 
September  1996  and  the  termination  of  it  due  to  a  power  loss  in  June  1997,  NSCAT  was 
very  successful  at  collecting  sea  surface  wind  vector  measurements.  Due  to  this  success, 
reliable  data  obtained  through  NSCAT  led  to  a  wide  variety  of  applications,  mostly 
addressing  scientific  and  operational  problems.  To  minimize  the  gap  between 
scatterometry  data,  the  schedule  of  the  following  QuikSCAT  mission  was  accelerated  and 
SeaWinds  was  launched  by  this  mission  in  June  1999 
(http://winds.jpl.nasa.gov/aboutScat/  history. cfm). 

2.  MISSION 

The  goal  of  the  QuikSCAT  mission  was  to  record  wind  speed  and  direction  data 
near  the  sea-surface  regardless  of  the  condition  of  the  weather  or  clouds.  It  was,  indeed,  a 
"quick  recovery"  mission  intended  to  reduce  the  gap  in  ocean-wind  vector  data  that 
resulted  from  the  loss  of  the  NSCAT  on  the  ADEOS.  QuikSCAT  was  launched  from  the 
Vandenberg  Air  Force  Base  in  California  onboard  of  a  vehicle  called  Titan  II.  The 
QuikSCAT  mission  was  conducted  for  NASA's  Earth  Science  Enterprise  by  a  division  of 
the  California  Institute  of  Technology  called  the  Jet  Propulsion  Laboratory  (JPL).  Other 
mission  partners  are  listed  in  Table  6. 
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WIND  SPEED  (m/s) 


Figure  25.  Hurricane  Dora  viewed  from  QuikSCAT  (From  JPL/NASA) 

The  orbit  that  QuikSCAT  was  designed  to  fly  in  was  a  circular,  near  polar  orbit 
that  was  limited  to  a  maximum  altitude  of  about  800  km  (500  miles)  above  Earth’s 
surface.  It  accomplished  a  full  orbit  in  around  101  minutes,  which  approximately 
corresponds  to  14  orbits  per  day.  SeaWinds,  which  is  an  active  radar  scatterometer,  is  the 
primary  instrument  on  the  QuikSCAT  satellite.  Each  day  SeaWinds  makes  around 
400,000  measurements  over  ocean,  land,  and  ice  continuously  with  a  1,800-kilometer¬ 
wide  band;  it  covers  90%  of  Earth’s  surface.  With  that  capability,  it  makes  hundreds  of 
times  more  observations  of  the  wind  velocity  on  the  Earth’s  surface  every  day  than  ships 
and  buoys  combined.  SeaWinds,  furthermore,  is  able  to  provide  measurements  of  both 
velocity  and  the  direction  of  the  surface  winds  that  are  high-resolution,  accurate, 
continuous,  and  which  do  not  depend  on  weather  conditions.  Obtaining  this  data  is 
crucial  in  terms  of  storm  warning,  global  climate  research,  and  many  more  different 
studies. 

There  were  several  ground  systems  involved  in  the  QuikSCAT  mission.  The  first 
types  were  Earth  Polar  Ground  Stations  that  tracked  the  satellite.  These  stations  were 
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located  in  Poker  Flats,  Alaska;  McMurdo,  Antarctica;  Svalbard,  Norway;  Wallops  Island, 
Virginia;  and  Hatoyama,  Japan  (the  back-up  station).  The  second  type  of  ground  system 
was  where  the  data  products  were  produced  and  distributed.  These  systems  were  the  Jet 
Propulsion  Laboratory  (JPL)  and  NOAA.  JPL  was  responsible  for  producing  high-quality 
research  data  products  and  distributing  them  to  the  science  community  within  two  weeks 
of  receipt.  This  distribution  was  made  through  the  JPL  Physical  Oceanography 
Distributed  Active  Archive  Center  (Po.daac).  NOAA,  on  the  other  hand,  was  responsible 
for  producing  operational  data  products  for  the  international  meteorological  community 
in  the  within  three  hours  of  receipt. 

Table  6.  Mission  Partners  of  QuikSCAT 

National  Oceanic  and  Atmospheric  Administration  (NOAA) 

NASA  Goddard  Space  Flight  Center  (GSFC) 

Ball  Aerospace  and  Technologies  Corporation 
U.S.  Air  Force  Space  and  Missile  Systems  Center 
Honeywell  Satellite  Systems  Operations 
Raytheon  E-Systems  Corporation 
Lockheed  Martin  Astronautics 
Hughes  Electron  Dynamics  Division 


QuikSCAT  was  in  service  until  there  was  a  mechanical  failure  of  the  spin 
mechanism  of  the  antenna.  Due  to  that  failure,  SeaWinds  stopped  rotating  in  November 
2009.  It  didn’t  stop  functioning  completely,  but  its  capability  was  significantly  reduced. 
However,  the  collected  data  have  still  been  setting  the  transfer  standard  for  cross¬ 
calibration,  and  the  standard  for  defining  the  measurement  stability  that  is  vital  for 
continuity  with  scatterometer  missions  in  the  future.  NASA  and  ISRO  have  agreed  to 
allow  the  QuikSCAT  team  to  have  access  to  measurements  taken  by  the  Oceansat-2 
scatterometer  (OSCAT),  which  has  been  operating  since  September  2009.  The 
partnership  between  NASA  and  ISRO  was  established  as  a  long-term  cooperation  that 
enables  direct  cross-calibration  of  QuikSCAT  with  OSCAT  to  help  in  the  generation  of 
an  ocean  vector  winds  climate  series  (http://winds.jpl.nasa.gov/missions/ 
quikscat/#mission). 
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IV.  METHODOLOGY 


A.  INTRODUCTION  TO  NEURAL  NETWORKS 

The  discovery  of  the  concept  of  neural  networks  is  a  captivating  progress  and  a 
very  important  milestone.  Frank  Rosenblatt  (1957)  was  the  inventor  of  the  first  neural 
network  for  pattern  classification.  This  network  is  called  “Perceptron.”  Thirteen  years 
after  its  invention,  Minsky  and  Papert  (1969)  published  Perceptrons,  which  called 
enormous  attention  to  neural  network  research.  The  concept  of  neural  networks  was 
formalized  and  demonstrated  in  this  book  as  well  as  discussing  a  few  of  its  serious 
limitations  concerning  the  original  architecture  (Noyes  1992).  To  be  more  specific, 
Perceptron,  in  this  publication,  is  proved  to  be  unable  to  execute  a  basic  logical 
computation  of  an  exclusive-or  (XOR).  That  disability  was  found  to  be  a  crucial 
drawback  at  first,  yet  research  continued,  and  as  a  result,  the  concept  of  neural  networks 
has  become  a  popular  tool  to  provide  solutions  to  complicated  problems  (Fitzpatrick 
1997;  Guthikonda  2005). 

The  neural  networks  approach  involves  modeling  the  biological  brain  with  the 
assistance  of  Artificial  Intelligence’s  brain-paradigm.  The  methodology  is  novel  as  it 
compares  the  previous  computer  science  concepts,  especially  in  terms  of  data  storage. 
Typical  programs  tend  to  store  the  data  in  solid  structure  types  and,  eventually,  the  data  is 
stored  in  a  centralized  database.  Neural  networks,  on  the  other  hand,  distribute  the 
information  throughout  the  network  (Noyes  1992;  Guthikonda  2005).  In  other  words, 
they  resemble  the  biological  brain. 

Neural  networks  attract  many  researchers  and  individuals  around  the  world  due  to 
their  promising  capability  to  solve  complex  problems.  They  can  be  extremely  fast  as  well 
as  efficient.  Such  features  enable  the  handling  of  large  data  sets  (Noyes  1992).  The 
reason  underlies  this  capability  is  that  each  node  contributing  to  a  neural  network  is 
fundamentally  self-governing.  Therefore,  only  a  small  portion  of  the  total  computation  in 
the  problem’s  grand-scheme  is  performed  by  each  node.  The  strength  lies  in  the 
aggregation  of  the  computational  power  of  all  the  nodes  in  the  entire  network.  This 
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special  design,  which  imitates  the  biological  brain,  makes  parallel  multitasking  feasible. 
Another  advantage  of  neural  networks  is  that  they  are  tolerable  against  fault  (Noyes 
1992).  This  means  that  a  limited  portion  of  corrupted  data,  or  a  limited  number  of  nodes 
that  exhibit  malfunction  or  do  not  even  function  at  all,  will  not  paralyze  the  whole 
network.  Under  these  circumstances,  the  network  figures  out  a  way  to  perform  reasonably 
even  with  these  faults  because  it  is  able  to  learn.  This  feature,  too,  is  inspired  by  the 
biological  brain  since  we  know  that  eevn  if  the  odrer  of  the  wrods  chnage,  we  are  sltil 
albe  to  usndranetd  a  stnecene  (this  sentence  with  misspelled  words  is  set  up  intentionally 
to  prove  that  the  statement  is  true)  (Guthikonda  2005). 

Before  delving  into  the  learning  capability  of  neural  networks,  the  definition  of 
learning  should  be  stated  clearly.  Learning  can  be  defined  in  two  ways  depending  on 
different  categories:  biologically  and  mechanically.  Biological  learning  is  detailed  as  an 
experience  which  causes  alteration  in  an  organism’s  state  and  improves  its  performance 
in  consecutive  similar  situations.  Mechanical  learning,  on  the  other  hand,  is  defined  as  a 
computational  approach  taken  in  order  to  achieve  new  intelligence  and  organize  that 
intelligence  in  order  to  acquire  new  abilities  (Noyes  1992).  In  light  of  these  two 
definitions,  it  can  be  assumed  that  a  neural  network  must  have  the  ability  to  learn  in  order 
to  be  considered  useful.  Training  of  neural  networks  leads  to  learning,  which  eventually 
enables  this  learning  to  be  applied  for  practical  purposes  (Guthikonda  2005). 

A  different  way  to  categorize  the  learning  process  is  to  take  control  on  the  data 
into  account.  Admittedly,  there  are  two  kinds  of  learning  in  this  regard:  supervised 
learning  and  unsupervised  learning  [Egge98].  In  supervised  learning,  the  accurate 
answers  are  known  and  the  main  purpose  is  to  train  the  network  to  respond  to  a  given 
situation  with  respect  to  these  accurate  answers.  Input  vectors  benefit  from  supervised 
learning  as  well  as  output  vectors.  Once  the  output  vectors  are  obtained,  they  are  used 
together  with  the  input  vectors  to  analyze  the  performance  of  the  network  and  to  locate 
any  errors.  Reinforcement  learning,  for  instance,  is  a  specialized  supervised  learning 
method,  which  only  provides  the  information  of  whether  or  not  the  output  is  accurate;  it 
is  based  on  the  back-propagation  algorithms  (Guthikonda  2005). 
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In  unsupervised  learning,  whereas,  it  is  either  the  accurate  answers  are  not 
available  in  advance  or  the  accurate  answers  are  not  shared  with  the  network  beforehand. 
The  input  vectors  are  handled  separately.  The  output  vectors  have  no  contribution  to  the 
learning  process,  unlike  in  supervised  learning.  The  network  must  perform  as  a  self- 
governed  system  and  must  extract  the  patterns  of  the  input  data  itself  without  an 
interaction  of  any  kind,  especially  human  interaction.  Supposedly,  this  is  a  crucial  feature 
since  it  is  time-consuming  and  sometimes  it  is  not  even  possible  for  a  human  to  make  the 
computations  by  him  or  herself  when  working  on  large  and  complicated  data  sets. 
Unsupervised  learning  is  the  kind  of  learning  utilized  by  self-organizing  maps 
(Guthikonda  2005). 

B.  THE  SELF-ORGANIZING  MAPS 

1.  INTRODUCTION 

Self-organizing  maps  were  first  described  by  Tuevo  Kohonen  in  1982.  “Self¬ 
organizing”  comes  from  the  fact  that  they  are  autonomous  and  supervision  is  not 
demanded.  SOMs  have  the  ability  to  learn  by  themselves  using  competitive  learning. 
“Maps”  refers  to  the  way  they  handle  the  given  input  data  sets,  which  is  by  attempting  to 
map  weights  of  the  data.  The  nodes  of  the  SOM  network  try  to  resemble  the  input  data 
sets  they  are  initially  given  and  this  reaction  of  the  nodes,  essentially,  sets  up  the 
foundation  for  the  whole  learning  process.  Preserving  the  main  and  important 
characteristics  of  the  input  data  sets  constitutes  the  fundamental  principle  of  SOMs,  and 
this  principle  is  what  makes  the  SOMs  special  in  comparison  to  other  methods.  Retaining 
the  topological  relationships  among  the  input  data  and,  mapping  these  relationships  to  an 
SOM  network  are  very  valuable  aspect  of  representing  complex  data  sets  (Guthikonda 
2005). 
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2. 


STRUCTURE 


SOMs  have  a  fairly  simple  structure,  which  is  illustrated  in  Figure  26. 


®  #  #  (§>  - 

Output  BMUs 


Figure  26.  Illustration  of  the  structure  of  a  SOM  (From  Maung  2012) 

A  4  x  4  SOM  network  is  shown  in  Figure  26  as  the  competition  layer.  There  are  a 
few  aspects  of  the  structure  that  are  worth  mentioning.  First  of  all,  each  input  node  has  a 
connection  with  each  map  node.  Even  for  a  small  network  (as  in  Figure  26)  that 
corresponds  to  a  large  number  of  connections.  For  example,  in  the  Figure  26  network,  4 
times  4  times  4  equals  64  connections.  Unlike  the  relationship  between  input  nodes  and 
map  nodes,  notice  that  there  is  no  connection  between  map  nodes  in  the  network.  The 
organization  of  the  nodes  is  a  two  dimensional  grid,  which,  as  mentioned  earlier,  provides 
an  easy  way  to  visualize  the  outputs.  This  type  of  structure  is  also  helpful  when  using  the 
SOM  algorithm.  Each  node  in  the  network  is  denoted  with  a  unique  coordinate  (i,  j).  This 
unique  coordinate  system  provides  an  easy  way  to  refer  to  a  node  in  the  network  and  to 
determine  the  ranges  between  nodes.  The  map  nodes  do  not  know  the  values  that  their 
neighbors  hold  due  to  the  fact  that  each  only  has  connections  with  the  input  nodes. 
Therefore,  each  node  must  use  the  information  given  by  the  input  vectors  in  order  to 
update  its  own  weights.  The  weight  vectors  of  both  the  map  nodes  and  the  input  vectors 

should  match  to  enable  the  algorithm  to  perform  properly  (Guthikonda  2005). 
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3. 


ALGORITHM 


As  an  unsupervised  learning  algorithm,  SOM  has  learning  and  prediction  phases. 
During  the  learning  phase,  map  is  built;  network  organizes  using  a  competitive  process 
with  a  training  set.  During  the  prediction  phase,  new  vectors  are  quickly  given  a  location 
on  the  converged  map,  easily  classifying  or  categorizing  the  new  data.  Vector 
quantization,  which  is  a  data  compression  technique,  is  the  other  crucial  property  of 
SOMs.  SOMs  reflect  multi-dimensional  data  sets  on  a  lower  dimensional  space;  in  most 
cases,  SOMs  reflect  these  data  sets  on  two  dimensions  (Buckland  2003).  This  provides  a 
user-friendly  and  easy-to-visualize  representation  of  the  complex  data  and  enables  a  user 
to  work  on  the  data  even  more  efficiently  (Guthikonda  2005).  The  algorithm  of  the  self¬ 
organizing  map  can  be  described  as  six  main  steps  (Buckland  2003).  These  steps  that  the 
algorithm  utilizes  are  outlined  in  Table  7  (Guthikonda  2005). 


Table  7.  Main  six  steps  of  the  algorithm  of  the  SOM 


1st  Step 

Each  node’s  weight  is  initialized  with  a  random  number  between  0  and  1. 

2nd 

Step 

A  random  vector  is  selected  from  the  training  data  set  and  introduced  to  the 
network. 

3rd  Step 

Each  node  in  the  network  is  inspected  in  order  to  find  out  which  one 
resembles  the  input  vector  more  accurately  in  terms  of  weights.  Best 

Matching  Unit  (BMU)  is  the  name  given  to  the  winning  node.  The  Euclidean 
distance  formula,  which  measures  the  similarity  between  two  data  sets,  is 
used  in  order  to  do  this  selection.  The  BMU  is  found  by  the  calculation  of  the 
distance  between  the  weights  of  node  and  the  input  vector,  (as  shown  in 
Equation  1). 

4th  Step 

The  radius  of  the  neighborhood  of  the  BMU  is  calculated.  This  radius  is 
initialized  as  the  network’s  radius  and  decays  with  each  time-step  until 
reaching  the  BMU  itself  (as  shown  in  Equations  2a  and  2b). 

5th  Step 

All  of  the  nodes  that  fall  in  the  radius  of  the  BMU  are  adjusted  in  order  to 
make  them  as  similar  as  possible  to  the  input  vector  (as  shown  in  Equations 

3a  and  3b).  The  weight  of  the  closest  node  to  the  BMU  is  changed  the  most, 
and  the  degree  of  this  change  diminishes  as  the  range  between  the  BMU  and 
the  nodes  increases  (as  shown  in  Equation  3c). 

6th  Step 

The  second  step  is  repeated  for  N  iterations. 
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Figure  27.  A  two-dimensional  illustration  that  shows  the  same  type  of  connection  of  each 
node  with  the  input  vectors.  Note  that  there  is  no  connection  among  the  nodes 
(Available  from  http://www.saedsayad.com/clustering_som.htm) 
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Figure  28.  An  illustration  that  shows  the  initialization  step  of  each  node’s  weights 
(From  http://www.saedsayad.com/clustering_som.htm) 
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The  observational  (input)  vector  V  =  (vi,  V2,  ...,  vn)  is  compared  to  the  node’s 
weight  vector  W  =  (w\,  vv’2,  . . wn).  The  distance  (or  difference)  between  the  two  vectors 
is  computed  by 


The  optimal  placement  of  the  weight  vectors  W (7)  is  obtained  from  the  minimization  of 
the  distance, 

D  =  min{Dy},  (2) 

j 

which  is  called  the  best  matching  unit  (BMU).  The  BMU  neighborhood  is  determined  as 
follows.  The  size  of  the  neighborhood  is  represented  by  an  exponential  decay  function 
that  shrinks  until  eventually  the  neighborhood  is  just  the  BMU  itself.  The  width  of  the 
lattice  at  iteration  t  (see  Figure  29)  a(t)  is  calculated  by 

^(0  =  cr0exP(—  -7),  (3) 

A 

where  00  is  the  width  of  the  initial  lattice;  X  is  a  constant.  The  new  weight  for  a  node  is 
the  old  weight,  plus  a  fraction  ( L )  of  the  difference  between  the  old  weight  and  the  input 
vector,  adjusted  (0)  based  on  distance  from  the  BMU, 

W(f  + 1)  -  W(f)  +  0(t)IXt)\Y(t)  -  W(f)] ,  (4) 

where  L(t)  is  the  learning  rate  at  the  iteration  t, 

L(t)  -  L0  exp(  ^-)  ,  (5) 

A 

with  Lo  the  initial  learning  rate;  and  0(t)  is  the  Gaussian  curve  for  the  neighborhood, 

r  d 2 

0(0=  exp - 5 —  ,  (6) 

2cr  (f)_ 

so  that  nodes  that  are  closer  are  influenced  more  than  farther  nodes.  As  iteration 
advances,  it  is  observed  that,  under  the  influence  of  this  formula,  a  given  node  weight 
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becomes  much  more  similar  to  the  input  vector  that  is  selected.  The  bigger  the  initial 
difference  between  a  node  and  selected  input  vector,  the  more  iteration  is  spent  to  train 
the  node  and  vice  versa.  That  difference  is  then  scaled  by  two  different  aspects:  the 
SOM’s  current  learning  rate  and  0(t)  (Guthikonda  2005). 
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Figure  29.  An  illustration  of  the  process  utilized  by  Equation  2a  (From  Guthikonda  2005) 


The  degree  of  influence  that  the  node’s  distance  from  the  BMU  has  on  the  node’s 
learning  is  shown  by  the  influence  rate.  Initially,  influence  rate  is  set  to  1  for  all  the 
neighboring  nodes  of  the  BMU  and  zero  for  the  ones  that  are  away  from  the  BMU.  At  the 
end,  due  to  weights’  random  distribution  and  numerous  iterations,  SOM  is  able  to  settle 
down  to  a  map  of  stable  zones  (http://www.saedsayad.com/clustering_som.htm).  Nodes 
that  do  not  fall  into  the  neighborhood  radius  are  entirely  ignored.  The  term 
“distFromBMU”  represents  the  exact  number  of  nodes  that  are  located  between  the  BMU 
and  the  given  node.  The  formula  used  to  calculate  that  particular  information  is: 

distFromBMU2  =  (bmul-nodel)2  +  (bmuJ-nodeJ)2,  (7) 
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The  fact  that  the  network  is  composed  of  nodes  that  are  on  a  two-dimensional  grid  makes 
the  calculation  possible.  The  amount  learned  by  the  nodes  on  the  edge  of  the 
neighborhood  radius  is  a  fractional  value  that  is  smaller  than  1.0.  Clearly  0(t)  has  a 
tendency  to  near  1 .0  if  the  term  “distFromBMU”  becomes  smaller.  In  terms  of  the  BMU, 
its  value  for  “distFromBMU”  will  be  zero;  this  enables  9(t)  to  reach  its  maximum  value, 
which  in  this  case  is  1.0  (Guthikonda  2005). 


Figure  30.  An  illustration  of  the  process  of  updating  BMU  along  with  its  neighbors 
towards  the  “x,”  which  in  this  case  is  the  input  sample  (From  Vesanto  et  al.  2000) 

The  equations  adopted  within  the  SOM  algorithm  vary  significantly  from  one 
study  to  another.  Besides  that,  the  biggest  debate  seems  to  be  regarding  the  optimal 
parameters.  Questions  concerning  the  sufficient  amount  of  iterations,  the  learning  rate, 
and  the  width  of  the  neighborhood  radius  have  still  not  been  answered  clearly.  The 
recommendation  of  Kohonen  in  this  matter  is  to  set  up  a  training  which  has  two  altered 
phases.  During  the  first  phase,  the  learning  coefficient  will  start  with  0.9  and  decrease  to 
0.1.  The  neighborhood  radius,  on  the  other  hand,  will  start  with  half  of  the  lattice’s 
diameter  and  decrease  to  the  closest  surrounding  nodes.  In  the  second  phase,  however, 
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the  learning  rate  will  start  with  0.1  and  decrease  to  0.0;  the  iteration  will  be  doubled  in 
comparison  with  the  first  phase.  The  value  for  the  neighborhood  radius  remains  as  1, 
which  encompasses  the  BMU  only.  In  the  light  of  these  descriptions,  Matthews  (2004) 
remarks  that  the  first  phase  can  be  denoted  as  “rough  tuning,”  which  enables  a  quick 
order;  and  the  second  phase  can  be  denoted  as  “fine  tuning,”  which  assures  a  precise 
representation  (Guthikonda  2005). 


Figure  31.  An  illustration  of  the  differences  among  the  four  neighborhood  functions. 
They  are,  from  left  to  right,  “bubble,”  “gaussian,”  “cutgass,”  and  “ep.”  The  first  row  is  a 
one-dimensional  representation,  whereas  the  second  row  is  two-dimensional  (From 

Vesanto  et  al.  2000). 
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Figure  32.  An  illustration  of  the  differences  among  the  three  functions  of  learning  rate. 
inv  is  represented  as  a  dashed  curve,  power  is  represented  as  a  dot-dashed  curve,  and 
linear  is  represented  as  solid  line  (From  Vesanto  et  al.  2000). 


4.  DAVIES  -  BOULDIN  INDEX 


The  Davies-Bouldin  index  (DBI) 


(8) 


is  exerted  to  the  data  sets  along  with  the  SOM,  to  determine  the  number  of  patterns  that 
best  represents  a  data  set.  DBI  was  introduced  by  David  L.  Davies  and  Donald  W. 
Bouldin  in  1979,  and  is  a  measure  that  is  used  for  assessing  clustering  algorithms.  DBI 
operates  as  an  internal  scheme,  and  the  authentication  of  how  precise  the  clustering  has 
been  accomplished  is  made  using  characteristics  of  the  data  set  (Davies  and  Bouldin 
1979).  In  the  formula;  N  is  the  number  of  clusters;  a;  is  the  average  distance  of  all 
patterns  in  cluster  i  to  their  cluster  centers  cp  oj  is  the  average  distance  of  all  patterns  in 
cluster  j  to  their  cluster  centers  cp  and  d(c;,  Cj)  is  the  distance  between  cluster  centers  c, 
and  Cj.  Small  values  of  DBI  correspond  to  clusters  that  are  compact,  and  whose  centers 
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are  far  away  from  each  other.  Consequently,  the  number  of  clusters  that  minimizes  DBI  is 
taken  as  the  optimal  number  of  clusters. 


5.  TOOLBOX  FOR  MATLAB 


1st  sample 
2nd  sample 
3rd  sample 
4th  sample 
etc. 
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Figure  33.  An  illustration  of  the  table  format  of  the  data  (From  Vesanto  et  al.  2000) 


The  Toolbox  can  perform  with  a  particular  form  of  data  which  is  called  table  data 
or  spreadsheet.  Table’s  each  row  represents  a  data  sample  and  is  composed  of  the 
variables  of  the  data  set.  These  variables  can  symbolize  different  aspects  of  the  data  set 
such  as  an  object’s  properties  or  measurements  taken  at  a  given  time.  The  crucial  rule 
governing  this  type  of  data  organization  is  that  samples  in  the  organization  must  share  the 
same  type  of  variables.  This  rule  assures  that  each  column  of  the  table  includes  all  values 
for  a  single  variable  (Vesanto  et  al.  2000). 

Either  symbolic  data  or  numeric  data  can  be  processed  with  the  Toolbox,  yet  the 
SOM  algorithm  only  allows  for  the  usage  of  numeric  data.  Furthermore,  the  order  of  the 
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objects  with  respect  to  the  values  they  contain  must  be  meaningful,  as  well.  This  means 
that  if  it  is  assumed  that  a  table  is  comprised  of  three  objects,  the  distance  between  the 
first  and  the  second  object  must  be  less  than  the  distance  between  the  second  and  third 
objects  (Vesanto  et  al.  2000). 

String  labels  of  the  data  samples  can  be  used  to  house  symbolic  data  in  the 
Toolbox,  which  enables  users  to  utilize  these  data  as  personal  notes  for  a  given  sample  in 
case  a  user  needs  to  remember  something  about  the  sample  or  s/he  needs  to  assure  a 
property  of  that  entity.  These  labels,  however,  will  not  be  recognized  by  the  SOM 
algorithm  as  mentioned  before.  Pyle  (1999)  proposes  that,  if  these  labels  are  aimed  to  be 
processed  in  the  training  algorithm,  conversion  methods  such  as  1-of-n  coding  or 
mapping  can  be  used  (Vesanto  et  al.  2000). 

C.  COMPARISON  WITH  OTHER  CONVENTIONAL  METHODS 

There  are  two  major  methods  often  preferred  when  identifying  patterns  of 
variability  in  oceanic  data  and  meteorological  data:  the  Empirical  Orthogonal  Function 
(EOF)  and  Principal  Component  Analysis  (PCA).  Liu  and  Weisberg  tested  both  EOF  and 
SOM  in  2005  and  2007  in  order  to  detect  patterns  of  ocean  currents  using  the  exact  same 
data  set.  Data  was  obtained  through  a  moored  Acoustic  Doppler  Current  Profiler  (ADCP) 
and  contained  a  long  velocity  time  series.  The  results  favored  SOM  as  researchers  found 
that  patterns  from  the  leading  mode  EOF  were  less  intuitive  and  less  accurate  compared 
to  the  SOM  patterns.  More  specifically,  non-symmetric  features  regarding  coastal  jet 
location,  current  strength,  and  veering  of  the  velocity  vector  with  respect  to  depth  that 
takes  place  between  the  current  patterns  of  upwelling  and  downwelling  were  detected  by 
the  nonlinear  SOM,  whereas  data  on  these  features  couldn’t  be  revealed  by  the  linear 
EOF  (Liu  and  Weisberg  2005;  Liu  and  Weisberg  2001). 

In  another  experiment,  SOM  and  EOF  were  compared  as  each  method’s 
extraction  of  an  artificial  data  that  represents  known  patterns  was  tested.  This  time  both 
SOM  and  EOF  were  found  to  be  successful  at  extracting  the  patterns  of  a  sine  wave  that 
is  linear  and  progressive,  but  SOM  was  also  successful  when  noise  was  added  to  the  data. 
On  the  other  hand,  when  the  experiment  was  carried  out  using  multiple  data  sets 
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concerning  more  complex  patterns,  the  EOF  technique  failed  to  detect  all  patterns  in 
contrast  with  the  SOM  method  which  was  successful  (Liu  et  al.  2006a;  Liu  and  Weisberg 
2001). 

Reusch  et  al.  (2005)  used  synthetic  data  sets  that  contained  both  negative  and 
positive  modes  of  four  sea  level  pressure  fields  of  the  North  Atlantic  that  are  idealized,  in 
order  to  test  PCA  method  against  SOM.  Test  is  executed  twice;  adding  and  disregarding 
noise  components.  At  the  end,  SOM  was  found  to  be  more  robust  in  comparison  to  PCA 
in  terms  of  the  extraction  of  the  patterns  of  variability  that  are  predefined.  Astel  et  al. 
(2007)  and  Annas  et  al.  (2007)  also  tested  PCA  method  against  SOM,  and  they  both 
confirmed  that  SOM  has  a  higher  performance  over  the  PCA  (Liu  and  Weisberg  2001). 

Another  artificial  neural  network  widely  used  and  preferred  in  terms  of  clustering 
is  K-means.  Bagao  et  al.  (2005)  compared  the  K-means  method  against  SOM,  and  came 
up  with  the  idea  that  SOMs  can  substitute  K-means  clustering  algorithms.  Lin  and  Chen 
(2006)  carried  out  another  test  to  compare  the  clustering  accuracy  among  the  SOM, 
Ward’s  method  and  the  K-means  method.  The  data  used  in  this  test  are  experimental  data 
sets  that  features  are  known  and  under  control  in  terms  of  cluster  membership  and  cluster 
dispersion.  SOM  method,  one  more  time,  found  to  be  more  accurate  in  specifying  the 
cluster  membership  compared  to  Ward’s  method  and  K-means  method  (Liu  and 
Weisberg  2001). 
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V.  SURFACE  CURRENTS 


The  geostrophic  velocity  data  from  1992  to  2012  were  obtained  from  AVISO  to 
provide  data  on  the  surface  currents.  Their  monthly  means  show  the  seasonal  variation  of 
the  surface  currents.  Six  patterns  were  found  after  conducting  the  SOM.  The  inter-annual 
variability  of  the  surface  currents  was  also  obtained.  Different  time  frames  were  selected 
and  examined  throughout  the  initial  analyses.  A  ten-year-long  data  set  was  found  to  be 
sufficient  when  extracting  the  patterns  of  the  long-term  temporal  and  spatial  variation  due 
to  the  fact  that  any  data  sets  covering  a  period  longer  than  10  years  do  not  affect  the 
results.  Therefore,  this  research  focuses  on  data  sets  from  1999  to  2009. 

A.  MONTHLY  MEAN  SURFACE  CURRENTS 

The  monthly  average  of  January  surface  currents  is  shown  in  Figure  34  As  shown 
in  the  figure,  many  eddies  with  various  sizes  are  evident.  The  Batumi  and  the  Sevastopol 
eddies,  both  of  which  are  anti-cyclonic,  are  the  most  pronounced  eddies  during  this 
month.  The  weak,  cyclonic  Suchumi  eddy  takes  place  north  of  the  Batumi  eddy.  During 
this  particular  time  of  the  year,  the  Caucasus  eddy  seems  immature,  but  there  is  a 
tendency  of  anti-cyclonic  circulation.  Off  of  the  Kerch  Strait,  the  well-organized,  anti- 
cyclonic  Kerch  eddy  is  detected.  Unlike  these  other  regions,  the  western  shelf  basin  is  not 
rich  in  terms  of  eddies.  Apart  from  the  Sevastopol  eddy,  the  anti-cyclonic  Bosphorus 
eddy  is  present.  Additionally,  the  northerly  western  boundary  current  is  distinctive.  Due 
to  the  abundance  of  eddies,  it  is  hard  to  identify  the  RIM  current,  but  the  general 
circulation  seems  to  have  two  main  gyres:  the  western  and  eastern  gyres. 

In  Figure  35,  the  climatological  average  of  the  surface  currents  for  February  is 
presented.  Most  of  eddies  that  existed  in  January  disappear  in  February.  The  anti-cyclonic 
Sevastopol  and  Batumi  eddies  still  exist  yet  they  are  weaker,  while  the  northerly  western 
boundary  current  is  not  detected.  The  main  western  and  eastern  gyres  are  observed  in 
terms  of  general  circulation. 
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Figure  35.  The  climatological  average  of  the  surface  currents  for  February 


In  March,  a  Batumi  dipole  eddies  is  detected  in  the  Batumi  eddy  region.  The 
Suchumi  eddy  appears  to  be  in  an  immature  stage,  but  there  is  a  tendency  of  cyclonic 
circulation.  The  Sevastopol  eddy  demonstrates  a  weak  cyclonic  circulation.  The  general 
circulation  of  the  basin  is  comprised  of  two  main  gyres,  with  the  eastern  gyre  being  more 
pronounced  (see  Figure  36). 

The  whole  basin  exhibits  chaotic  circulation  patterns  in  April.  The  Sevastopol 
eddy  disappears  and  the  same  Batumi  dipole  eddies  is  evident  in  the  Batumi  eddy  region. 
The  western  main  gyre  is  anti-cyclonic  whereas  the  eastern  main  gyre  is  cyclonic.  Both 
main  gyres  are  strong  during  this  period  (see  Figure  37). 

In  May,  the  Batumi  dipole  eddies  in  the  southeastern  basin  starts  detaching  from 
the  coast  and  gets  weaker,  as  does  the  anti-cyclonic  Caucasus  eddy.  A  strong,  northerly 
boundary  layer  arises  in  the  southeastern  part  of  the  basin.  During  this  time  of  the  year, 
both  the  western  and  eastern  main  gyres  are  cyclonic  and  well-defined  (see  Figure  38). 

In  June,  the  whole  basin  is  dominated  by  eddies  with  various  sizes  and  strength. 
The  Batumi  dipole  eddies  in  the  Batumi  eddy  region  detaches  completely  from  the  coast 
and  moves  to  the  north.  Two  small  cyclonic  eddies  are  evident  in  the  Sevastopol  eddy 
region  and  the  Sakarya  eddy  begins  to  develop  (see  Figure  39). 

Two  cyclonic  main  gyres  are  evident  during  July.  The  Bosphorus  eddy  is  detected 
at  its  developing  stage.  The  Batumi  dipole  eddies  in  the  southeastern  part  of  the  basin 
disappears  and  the  Batumi  eddy  appears  as  a  single  eddy  structure.  The  Sevastopol  eddy 
seems  to  be  included  into  the  Western  main  gyre.  The  southerly  and  northerly  boundary 
layers  are  formed  in  the  western  and  eastern  parts,  respectively  (see  Figure  40). 

The  general  circulation  is  formed  by  an  elongated  cyclonic  single  loop  called  the 
RIM  current  in  August.  Although  there  are  regions  showing  a  tendency  of  a  circular 
motion,  the  only  well-pronounced  eddy  is  the  cyclonic  Batumi  eddy.  The  anti-cyclonic 
tendency  next  to  the  Batumi  eddy  is  considered  a  sign  of  the  birth  of  the  Suchumi  eddy 
(see  Figure  41). 

The  cyclonic  RIM  current  persists  in  September,  as  well.  During  this  time  of  the 
year  the  Batumi  eddy  seems  to  get  stronger  and  wider.  Meanwhile,  the  Sevastopol  eddy  is 

67 


still  not  evident  and  continues  to  be  part  of  the  general  circulation.  Neither  major  nor 
minor  eddies  are  detected  in  the  open  basin  (see  Figure  42). 

In  October,  the  Batumi  eddy  gets  stronger  and  moves  to  the  southeastern  corner  of 
the  basin.  The  elongated  single  loop  structure  appears  to  be  dissolved  and  two-main-gyre 
structure  starts  to  take  over  again.  During  this  time,  the  eastern  main  gyre  is  stronger  than 
the  western  main  gyre.  The  Kerch  and  the  Caucasus  eddies  start  to  develop  in  the 
northeastern  part  of  the  basin  (see  Figure  43). 


Longitude  (°E) 

Figure  36.  The  climatological  average  of  the  surface  currents  for  March 
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Figure  37.  The  climatological  average  of  the  surface  currents  for  April 


Figure  38.  The  climatological  average  of  the  surface  currents  for  May 
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Figure  39.  The  climatological  average  of  the  surface  currents  for  June 


Figure  40.  The  climatological  average  of  the  surface  currents  for  July 
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Figure  42.  The  climatological  average  of  the  surface  currents  for  September 


Figure  44.  The  climatological  average  of  the  surface  currents  for  November 
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Figure  45.  The  climatological  average  of  the  surface  currents  for  December 

During  the  month  of  November,  the  Batumi  eddy  increases  its  strength,  affecting 
a  larger  area  in  the  southeastern  corner  of  the  basin.  Meanwhile,  the  two-main-gyre 
structure  persists.  Both  main  gyres  are  strong  and  cyclonic  during  this  time  of  the  year.  In 
November  the  Sevastopol  eddy  disappears  completely  and,  in  the  southern  region,  the 
Kizilirmak  eddy  starts  to  develop  (see  Figure  44). 

In  December,  the  anti-cyclonic  Batumi  eddy  is  at  its  maximum  strength.  In  the 
northwestern  region  the  Sevastopol  eddy  forms  again  and  rotates  anti-cyclonically.  The 
main  gyres  are  very  weak  on  both  sides  of  the  basin.  In  the  northeastern  and  eastern 
regions,  the  Kerch  eddy  and  the  Suchumi  eddy,  respectively,  develop  (see  Figure  45). 
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B.  TEMPORAL  AND  SPATIAL  VARIABILITY 


To  study  the  temporal  and  spatial  variation  of  the  surface  currents,  the  Black  Sea 
basin  is  subject  to  iterative  SOM  analysis.  The  DBI  demonstrates  the  best  option  with 
regards  to  the  selection  of  total  number  of  clusters  to  represent  the  whole  data  set.  The 
smaller  the  DBI,  the  better  and  more  practical  the  data  representation  is.  Therefore,  in 
light  of  the  DBI,  the  SOM  analysis  for  the  surface  currents  is  based  on  six  clusters  (see 
Figure  46). 
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Figure  46.  The  Davies-Bouldin  index  results  concerning  the  geostrophic  velocity  data  set 

1.  SIX  PATTERNS 

a.  Pattern  1:  Sevastopol  Cyclonic  and  Batumi  Dipole  Eddies 

The  first  pattern  represents  more  than  20%  of  the  whole  data  set.  In  this 
pattern,  the  Batumi  dipole  eddies  are  detected  at  the  southeastern  corner  of  the  basin.  The 
eddy  on  the  right  is  strong,  cyclonic,  and  it  spins  at  ~8  cm/s,  whereas  the  other  eddy  is 
weak  and  anti-cyclonic.  At  the  northwestern  comer  of  the  basin,  the  cyclonic  Sevastopol 
eddy  which  spins  at  ~5  cm/s  is  evident.  The  general  circulation  is  comprised  of  two  main 
gyres.  The  western  main  gyre  is  anti-cyclonic  and  weak.  The  eastern  main  gyre,  on  the 
other  hand,  is  cyclonic  and  strong.  The  northerly  boundary  current  forms  in  the  west.  No 
eddy  formation  is  observed  in  the  open  parts  of  the  basin  (see  Figure  47). 
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Geostrophic  Velocity  Pattern  :  1  (20.8333%) 
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The  first  pattern  of  the  surface  currents  and  its  percentage 


Geostrophic  Velocity  Pattern  :  2  (16.6667%) 
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Figure  48.  The  second  pattern  of  the  surface  currents  and  its  percentage 


b.  Pattern  2:  Cyclonic  RIM  Current  and  Anti-cyclonic  Batumi  Eddy 

More  than  16%  of  the  data  set  is  represented  by  the  second  pattern.  The 
main  feature  of  this  pattern  is  the  strong,  anti-cyclonic  Batumi  eddy  spinning  at  ~8  cm/s. 
Apart  from  the  Batumi  eddy,  no  other  eddy  structure,  including  the  major  Sevastopol 
eddy,  is  observed.  The  general  circulation  is  formed  by  the  strong  and  cyclonic  RIM 
current  flowing  at  ~5  cm/s.  The  open  parts  of  the  basin  are  relatively  calm  as  opposed  to 
the  coastal  regions.  Therefore,  this  particular  pattern  is  mostly  dominated  by  the  RIM 
current  (see  Figure  48). 

c.  Pattern  3:  Anti-cyclonic  Sevastopol  and  Batumi  Eddies 

The  third  pattern  also  represents  almost  16%  of  the  whole  data  set.  The 
major  eddies  of  the  basin  are  the  Batumi  and  Sevastopol  eddies  which  spin  at  ~7-8  cm/s. 
Both  of  these  major  eddies  are  strong  and  anti-cyclonic.  The  RIM  current  does  not  appear 
which  means  that  the  general  circulation  is  formed  by  the  western  and  eastern  main 
gyres.  In  the  northeastern  corner,  the  weak  Caucasus  eddy  exists,  while  the  strong 
northerly  boundary  current  is  present  in  the  west  (see  Figure  49). 

d.  Pattern  4:  Cyclonic  RIM  Current  and  Batumi  Eddy 

The  fourth  pattern  comprises  more  than  20%  of  the  whole  representation 
of  the  data  set.  In  this  pattern  the  Black  Sea  basin  is  almost  entirely  dominated  by  the 
strong  cyclonic  RIM  current  which  flows  at  ~10  cm/s.  The  open  parts  of  the  sea  stay 
relatively  calm  and  the  very  weak  cyclonic  Batumi  eddy  is  detected  but  it  is  nearly 
absorbed  by  the  RIM  current  (see  Figure  50). 

c.  Pattern  5:  Anti-cyclonic  RIM  Current  and  Batumi  Dipole  Eddies 

The  fifth  pattern  represents  15%  of  the  whole  data  set  and  it  demonstrates 
a  basin- wide,  chaotic  environment.  At  the  southeastern  corner  the  Batumi  dipole  eddies 
take  place  and  the  stronger  eddy  spins  at  ~7  cm/s.  The  weak,  anti-cyclonic  Crimea  eddy 
is  observed  in  the  north,  whereas  the  Sevastopol  eddy  disappears.  The  general  circulation 
is  formed  by  the  strong,  anti-cyclonic  RIM  current  which  flows  at  -6-1  cm/s  (see  Figure 
51). 
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Geostrophic  Velocity  Pattern  :  3  (15.8333%) 
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Figure  50.  The  fourth  pattern  of  the  surface  currents  and  its  percentage 
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Figure  51.  The  fifth  pattern  of  the  surface  currents  and  its  percentage 

a.  Pattern  6:  Anti-cyclonic  RIM  Current  and  Multi  Eddies 

More  than  10%  of  the  total  representation  for  the  data  set  is  covered  by  the 
sixth  pattern.  The  southeastern  comer  of  the  basin  is  entirely  dominated  by  the  strong, 
anti-cyclonic  Batumi  eddy  spinning  at  -8  cm/s.  At  the  north  of  the  Batumi  eddy,  the 
cyclonic  Suchumi  eddy  is  detected.  The  weak,  anti-cyclonic  Kerch  eddy  is  observed  in 
the  north.  At  the  northwestern  corner  the  strong  anti-cyclonic  Sevastopol  eddy  forms, 
which  spins  at  ~6  cm/s.  The  general  circulation  is  comprised  of  the  strong,  anti-cyclonic 
RIM  current  (see  Figure  52). 


Geostrophic  Velocity  Pattern  :  5  (15%) 
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Figure  52.  The  sixth  pattern  of  the  surface  currents  and  its  percentage 

2.  VORTICITY  FEATURES 

The  33.7°  E  longitude  is  taken  to  divide  the  Black  Sea  into  eastern  and  western 
parts.  The  spatially  averaged  Pattern-/  relative  vorticity  (Q  and  enstrophy  (is,-)  for  the  two 
parts  are  calculated: 

Ce  =^jj  Cfdxdy,  Cw  =^-JJ  C.dxdy,  (9) 

4" = 4- If  £dxdy-  4' = ^  If  ft no) 

where  the  subscript  ‘£”  means  the  East  Black  Sea,  and  ‘IT’  indicates  the  West  Black  Sea. 

The  enstrophy,  basically,  describes  how  the  fluid  acts  to  dissipate  energy  in  terms  of  The 

The  enstrophy,  basically,  describes  how  the  fluid  acts  to  dissipate  energy  in  terms  of 

kinetic  energy.  In  a  mechanical  problem,  this  would  be  like  how  friction  acts  to  slow 

down  movement  (velocity)  because  it  dissipates  energy.  The  averaged  enstrophies  over 
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Geostrophic  Velocity  Pattern  :  6  (10.8333%) 
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the  six  patterns  for  the  east  and  West  Black  Sea  were  calculated,  and  the  differences  of 
each  pattern  versus  the  average  were  taken  as  the  anomalies  relative  to  the  mean 
enstrophies.  Figure  53  shows  the  scatter  diagrams  of  each  pattern  in  terms  of  the  East  and 
West  Black  Sea  relative  vorticity  and  enstrophy  anomaly. 
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Figure  53.  Horizontally  averaged  East  and  West  (a)  relative  vorticity,  and  (b)  enstrophy 
anomaly  values  of  the  six  patterns  of  the  surface  currents 


3.  TEMPORAL  VARIATION 

To  better  illustrate  the  monthly  and  seasonal  variability  of  the  surface  currents, 
the  rate  of  occurrence  of  the  six  patterns  was  computed  and  the  results  are  shown  in 
Figure  54.  Pattern  1  is  one  of  the  two  patterns  that  dominate  during  winter-like  and  fall¬ 
like  months.  Its  maximum  contribution  is  in  January  with  -45%  and  it  disappears 
completely  in  June.  Even  though  there  is  less  than  a  35%  chance  it  will  appear,  it  is  the 
pattern  most  likely  to  be  observed  in  December.  Pattern  2  is  more  likely  to  be  observed 
during  the  spring.  Its  maximum  contribution  is  detected  in  April  with  more  than  50%  and 
it  has  the  highest  percentage  among  the  six  patterns  in  this  particular  month.  In  summer 
and  fall  it  tends  to  stay  at  around  20%,  it  declines  in  winter-like  months,  and  it 
completely  disappears  in  June.  Pattern  3  tends  to  appear  towards  the  end  of  the  spring  and 
the  beginning  of  the  summer.  It  reaches  its  maximum  in  May  and  becomes  the  pattern 
most  likely  to  appear  during  that  time  of  the  year.  It  demonstrates  an  oscillation  during 
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fall-like  months  and  never  shows  up  in  January  and  February.  The  first  three  months  of 
the  year  are  dominated  by  the  Pattern  4.  It  tends  to  appear  strongly  during  cold  periods, 
whereas  during  warm  periods  it  exhibits  a  weak  contribution.  Pattern  4  disappears  in  May 
and  doesn’t  appear  again  until  September.  Despite  the  fact  that  Pattern  5  shows  maximum 
percentages  in  fall-like  months  and  dominates  in  October,  it  has  no  significant  seasonality 
and  occurs  occasionally  throughout  the  year.  Pattern  6  completely  disappears  from 
November  to  April,  and  appears  only  five  months  during  the  year.  It  dominates  the  period 
from  June  to  September  and  reaches  its  maximum  in  July.  With  more  than  40% 
contribution,  it  becomes  the  pattern  most  likely  to  be  detected  in  July. 
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Figure  54.  Monthly  percentage  of  the  six  patterns  of  the  surface  currents 


The  time  series  of  the  six  patterns  of  the  surface  currents  from  2000  through  2009, 
with  one  year  overlap,  that  was  derived  by  SOM  is  shown  in  Figure  55  and  the  evolution 
of  the  duration  and  the  frequency  of  the  same  patterns  is  shown  in  Figure  56  in  order  to 
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illustrate  the  inter-annual  variability.  Furthermore,  Figure  57,  Figure  58,  and  Figure  59 
demonstrate  the  contemporaneous  EAWR,  NAO,  and  ENSO  indices  respectively,  along 
with  the  times  series  of  the  surface  currents.  In  addition  to  those  figures,  Table  9  is 
provided  to  allow  a  deeper  understanding  of  the  relationship  between  the  surface  currents 
and  the  large-scale  teleconnection  indices. 

The  potential  influence  of  large-scale  teleconnection  indices  on  the  recurrent 
patterns  of  the  surface  currents  (shown  in  Figures  47-52)  are  evaluated  by  estimating  for 
each  recurrent  pattern  the  value  of  the  climate  indices  (i.e.  EAWR,  NAO,  ENSO), 
acquired  for  the  different  months  mapped  towards  this  pattern.  This  process  enables 
determining  whether  a  recurrent  pattern  is  associated  to  either  positive  or  negative  value 
of  a  particular  index  significantly,  and  therefore  a  connection  between  the  surface  winds 
over  the  Black  Sea  and  large-scale  teleconnection  indices.  Figure  60  shows  for  each 
recurrent  pattern  of  the  surface  current,  the  distribution  of  positive  and  negative  phases  of 
the  three  teleconnection  indices  during  the  months  mapped  towards  this  atmospheric 
pattern.  According  to  the  figure,  the  EAWR  index  and  the  ENSO  index  are  negative 
significantly  during  months  mapped  towards  pattern  5,  which  indicates  that  the 
occurrence  of  pattern  5  (Anti-cyclonic  RIM  Current  and  Batumi  Dipole  Eddies)  is 
promoted  by  these  indices.  On  the  other  hand,  ENSO  index  is  significantly  positive 
during  months  mapped  towards  pattern  2  (Cyclonic  RIM  Current  and  Anti-cyclonic 
Batumi  Eddy)  and  therefore  influence  the  surface  current  structure. 
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Geostrophic  Velocity  Patterns  from  2000  to  2005 


Geostrophic  Velocity  Patterns  from  2004  to  2009 


Figure  55.  The  inter-annual  variability  of  the  six  patterns  of  the  surface  currents 


Year 

Figure  56.  Evolution  of  the  duration  and  the  frequency  of  the  six  patterns  of  the  surface 
currents  shown  in  Figures  47  -  52.  For  instance,  the  SOM  Pattern  2  lasted  for 
approximately  3  months  and  dominated  the  period  between  2002  and  2003  when  no  other 

patterns  occurred. 
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Geostrophic  Velocity  Patterns  from  2000  to  2009 


Figure  57.  Comparison  between  the  inter- annual  variability  of  the  six  patterns  of  the 

surface  currents  and  the  EAWR  index 


Geostrophic  Velocity  Patterns  from  2000  to  2009 


Figure  58.  Comparison  between  the  inter- annual  variability  of  the  six  patterns  of  the 

surface  currents  and  the  NAO  index 
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Geostrophic  Velocity  Patterns  from  2000  to  2009 


Figure  59.  Comparison  between  the  inter-annual  variability  of  the  six  patterns  of  the 

surface  currents  and  the  ENSO  index 


Table  8.  The  percentages  of  the  six  patterns  of  the  surface  currents  with  respect  to 
the  large-scale  teleconnection  indices.  For  example,  when  the  ENSO  index 
is  positive,  Pattern  2  has  the  highest  possibility  to  form  among  the  six 
patterns  at  25%.  Yet,  when  the  ENSO  index  is  negative,  pattern  6  has  the 
smallest  possibility  of  occurring  at  only  8%. 


GEOSTROPHIC  VELOCITY 

(+)  NAO 

(-)  NAO 

(+)  EAWR 

(-)EAWR 

(+)ENSO 

(-)  ENSO 

Pattern  1 

20% 

13  % 

20% 

14% 

11  % 

20% 

Pattern  2 

25  % 

15  % 

25  % 

16  % 

30% 

11  % 

Pattern  3 

14% 

16  % 

11  % 

18  % 

14.5  % 

14.5  % 

Pattern  4 

15  % 

24% 

22  % 

16  % 

20% 

23  % 

Pattern  5 

17  % 

17  % 

12  % 

22  % 

7  % 

23.5  % 

Pattern  6 

9% 

15  % 

10% 
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17.5  % 

8  % 
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Figure  60.  Distribution  for  each  pattern  of  the  surface  currents  between  positive  (red)  and 
negative  (blue)  phases  of  the  teleconnection  indices.  The  areas  with  darker  colors 
represent  the  strong  phases  for  which  the  index’s  absolute  value  is  greater  than  a  standard 
deviation.  Monthly  average  values  of  the  climate  are  obtained  through  the  NOAA 
Climate  Prediction  Center  (Available  at  http://www.cpc.ncep.noaa.gov/data/teledoc/ 

telecontents .  shtml) . 
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VI.  SURFACE  WINDS 


The  wind  data  from  1992  to  2012  were  obtained  from  QuikSCAT.  This  data  is 
representative  of  the  winds  near  the  surface  of  the  Black  Sea.  The  data  provides  monthly 
means  that  show  the  seasonal  variation  of  the  surface  currents.  After  conducting  the 
SOM,  six  patterns  were  found.  Additionally,  the  inter-annual  variability  of  the  winds  was 
obtained.  Throughout  the  initial  analysis  different  time  frames  were  selected  and 
examined  throughout  the  initial  analyses.  A  ten-year-long  data  set  was  found  to  be 
sufficient  to  extract  the  patterns  of  the  long-term  temporal  and  spatial  variation  due  to  the 
fact  that  any  data  sets  covering  a  period  longer  than  ten  years  do  not  affect  the  results. 
Therefore,  this  research  focuses  on  data  sets  from  1999  to  2009. 

A.  MONTHLY  MEAN  WINDS 

The  monthly  average  of  the  winds  for  January  is  shown  in  Figure  61.  In  the  Sea  of 
Azov  region,  the  winds  are  weak  and  easterly  with  a  slight  cyclonic  curve.  Strong 
northerly  winds  are  observed  in  the  western  shelf  and  they  get  weaker  as  they  move 
towards  the  central  Black  Sea.  At  the  central  eastern  shelf,  a  weak  cyclonic  motion  of  the 
winds  is  detected.  These  winds  are  very  weak  compared  to  the  winds  of  the  western  shelf. 
Meanwhile,  the  wind  in  the  Batumi  eddy  region  shows  a  strong  southerly  pattern  within  a 
small  region.  At  the  southeastern  corner  of  the  basin,  weak  easterly  winds  are  observed. 

In  Figure  62,  the  climatological  average  of  the  winds  for  February  is  presented. 
This  average  is  similar  to  January  with  some  exceptions.  Additionally,  the  winds  in 
February  are  stronger.  During  this  time,  the  Sea  of  Azov  exhibits  a  chaotic  pattern 
comprised  of  very  weak  winds,  the  cyclonic  motion  at  the  eastern  shelf  is  more 
pronounced,  the  southerly  winds  in  the  Batumi  eddy  region  become  a  part  of  the  cyclonic 
motion,  and  the  weak  easterly  winds  at  the  southeastern  comer  of  the  basin  disappear. 

In  March,  the  winds  in  the  Sea  of  Azov  region  become  stronger  as  they  begin  to 
blow  from  the  southeast.  The  winds  north  of  the  western  basin  take  on  a  chaotic  nature  as 
their  movement  appears  to  have  no  apparent  pattern.  The  winds  at  the  western  and 
eastern  shelves  are  weaker.  The  central  Black  Sea  region  exhibits  relatively  strong  winds. 
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The  cyclonic  motion  of  the  winds  at  the  eastern  shelf  remains  and  the  strong  winds  in  the 
Batumi  eddy  region  appear  again  (see  Figure  63). 


Figure  64  shows  the  climatological  average  of  the  winds  for  April.  A  clear  pattern 
can  be  observed  as  the  strong  northerly  winds  start  off  of  the  Crimean  peninsula  and  blow 
through  the  central  basin  reaching  to  the  Bosphorus.  No  clear  pattern  is  observed  at  the 
eastern  shelf  except  for  the  mediocre  southerly  winds  at  the  southeast  corner.  Weak 
easterly  winds  are  observed  at  the  Sea  of  Azov. 


Figure  61.  The  climatological  average  of  the  winds  for  January 
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Figure  62.  The  climatological  average  of  the  winds  for  February 
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Figure  63.  The  climatological  average  of  the  winds  for  March 


Figure  64.  The  climatological  average  of  the  winds  for  April 


In  May  the  strong  northerly  winds  are  evident  at  the  north  of  the  eastern  shelf  and 
at  the  south  of  the  western  shelf.  They  curve  slightly  anti-cyclonically  at  the  southwest  of 
the  western  shelf.  Whereas,  the  winds  at  southeast  of  the  eastern  shelf  exhibit  a  cyclonic 
curve.  A  very  weak  anti-cyclonic  circulation  can  be  noted  at  the  northwestern  corner  of 
the  basin.  Additionally,  weak  northerly  winds  are  detected  at  the  Sea  of  Azov 
(see  Figure  65). 

Figure  66  showcases  the  climatological  average  of  the  winds  for  June.  The  winds 
mostly  have  a  basin-wide  northwesterly  pattern  and  are  very  strong.  A  slight  anti- 
cyclonic  curve  at  the  southwestern  comer  of  the  western  shelf  is  evident.  At  the  far 
southeast,  a  tendency  of  westerly  winds  is  noted. 

Figure  67,  Figure  68,  Figure  69,  and  Figure  70  illustrate  the  climatological 
averages  of  the  winds  for  July,  August,  September,  and  October,  respectively.  The  nature 
of  the  wind  patterns  is  very  similar  for  the  climatological  averages  of  these  months  with 
minor  exceptions.  The  winds,  in  general,  become  northerlies  and  get  stronger  basin-wide 
in  July.  During  these  months,  the  anti-cyclonic  curve  in  the  western  shelf  and  the 
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cyclonic  curve  in  the  eastern  shelf  are  well-pronounced.  In  August  the  winds  get  even 
stronger  and  they  differentiate  from  the  winds  in  July  in  that  the  winds  at  the  Sea  of  Azov 
start  to  blow  from  the  northeast.  The  wind  pattern  of  August  is  the  same  in  September 
except  that  the  winds  weaken  at  the  Sea  of  Azov.  In  October  the  winds  generally  turn 
slightly  to  the  west  from  their  northerly  flow  pattern. 

In  Figure  71  the  climatological  average  of  the  winds  for  November  is  presented. 
The  winds  start  to  flow  from  the  northeast  instead  of  from  the  north  and  they  get  weaker 
at  the  Sea  of  Azov.  The  northerly  winds  dominate  the  southern  part  of  the  western  shelf, 
the  northern  part  of  the  eastern  shelf,  and  the  central  Black  Sea.  Meanwhile,  a  weak 
cyclonic  motion  at  the  far  southeast  is  detected. 

In  December  only  the  winds  in  the  southern  part  of  the  eastern  shelf  remain  weak, 
while  the  winds  in  the  rest  of  the  basin  get  stronger.  The  northerly  pattern  of  the  winds 
continues  to  dominate.  Additionally,  northeasterly  winds  are  noted  at  the  Sea  of  Azov 
and  the  motion  of  the  winds  at  the  southeastern  comer  of  the  basin  appears  to  be  chaotic 
rather  than  cyclonic  (see  Figure  72). 


Figure  65.  The  climatological  average  of  the  winds  for  May 
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Figure  66.  The  climatological  average  of  the  winds  for  June 
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Figure  67.  The  climatological  average  of  the  winds  for  July 
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Figure  68.  The  climatological  average  of  the  winds  for  August 
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Figure  69.  The  climatological  average  of  the  winds  for  September 
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Figure  70.  The  climatological  average  of  the  winds  for  October 
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Figure  71.  The  climatological  average  of  the  winds  for  November 
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Figure  72.  The  climatological  average  of  the  winds  for  December 


B.  TEMPORAL  AND  SPATIAL  VARIABILITY 

To  study  the  temporal  and  spatial  variation  of  the  winds,  the  Black  Sea  basin  is 
subject  to  iterative  SOM  analysis.  The  DBI  demonstrates  the  best  option  in  terms  of 
representing  the  total  number  of  clusters  in  the  data  set.  The  smaller  the  DBI,  the  better 
and  more  practical  the  data  representation  is.  Therefore,  in  light  of  the  DBI,  the  SOM 
analysis  for  the  surface  currents  is  based  on  six  clusters  (see  Figure  73). 
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Figure  73.  The  Davies-Bouldin  index  results  concerning  the  wind  data  set 
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1. 


SIX  PATTERNS 


a.  Pattern  1:  The  Northerly  Winds  and  the  Cyclonic  Curve  in  the 
Western  Shelf 

The  first  pattern  represents  more  than  13%  of  the  whole  data  set.  In  this 
pattern,  westerly  winds  are  observed  north  of  the  western  shelf.  The  winds  curve  anti- 
cyclonically  and  become  northeasterly  south  of  the  western  shelf.  At  the  eastern  self,  on 
the  other  hand,  the  winds  blow  from  the  northwest  with  a  slight  cyclonic  curve  south  of 
the  eastern  shelf  (see  Figure  74). 

b.  Pattern  2:  The  Easterly  Winds  and  the  Anti-cyclonic  Curve  in  the 
Western  Shelf 

More  than  13%  of  the  total  representation  for  the  data  set  is  covered  by  the 
second  pattern.  The  main  feature  of  this  pattern  is  the  anti-cyclonic  curve  of  the  winds  in 
the  western  shelf.  The  winds  blow  from  the  east,  starting  at  the  eastern  shelf,  and  from 
the  central  Black  Sea  region.  The  winds  are  stronger  to  the  west  and  they  are  strongest 
north  of  the  western  shelf  (see  Figure  75). 

c.  Pattern  3:  The  Dominance  of  the  Strong  Northerly  Winds 

The  third  pattern  represents  23%  of  the  whole  data  set.  The  strong 
northerly  winds  almost  completely  dominate  the  Black  Sea  region.  The  winds  lose  their 
strength  to  the  south  of  the  eastern  shelf  and  they  exhibit  a  cyclonic  curve  in  this  region 
(see  Figure  76). 
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Figure  74.  The  first  pattern  of  the  winds  and  its  percentage 


Geostrophic  Wind  Pattern  :  2  (13.6752%) 
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Figure  75.  The  second  pattern  of  the  winds  and  its  percentage 
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d.  Pattern  4:  The  Northwesterly  Winds  in  the  West  and  the  Westerly 

Winds  in  the  East 

The  fourth  pattern  comprises  more  than  16%  of  the  whole  representation 
of  the  data  set.  The  western  and  eastern  shelves  exhibit  two  different  wind  characteristics. 
The  strong  northwesterly  winds  occur  in  the  western  shelf  with  a  little  anti-cyclonic  curve 
to  the  south,  whereas,  starting  from  the  central  Black  Sea,  the  winds  blow  from  the  west 
and  the  whole  eastern  shelf  is  dominated  by  the  relatively  weak  easterly  winds 
(see  Figure  77). 

e.  Pattern  5:  The  Dominance  of  the  Strong  Northeasterly  Winds 

The  fifth  pattern  represents  more  than  21%  of  the  whole  data  set.  The 
northeasterly  winds  dominate  the  whole  Black  Sea  region  with  the  exception  of  the 
southeastern  comer  of  the  basin.  The  winds  become  weaker  to  the  southeast  and  almost 
disappear  completely  close  to  the  coast  (see  Figure  78). 


Geostrophic  Wind  Pattern  :  3  (23.0769%) 
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Figure  76.  The  third  pattern  of  the  winds  and  its  percentage 
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Geostrophic  Wind  Pattern  :  4  (16.2393%) 


28  30  32  34  36  38  40 


Longitude  (°E) 

Figure  77.  The  fourth  pattern  of  the  winds  and  its  percentage 


Geostrophic  Wind  Pattern  :  5  (21.3675%) 
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Figure  78.  The  fifth  pattern  of  the  winds  and  its  percentage 
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Figure  79.  The  sixth  pattern  of  the  winds  and  its  percentage 

/.  Pattern  6:  The  Strong  Southwesterly  Winds 

Almost  12%  of  the  total  representation  for  the  data  set  is  covered  by  the 
sixth  pattern.  The  southwesterly  winds  dominate  the  whole  Black  Sea  region  with  the 
exception  of  the  southeastern  corner  of  the  basin.  The  winds  blowing  from  the  southeast 
become  weaker  and  they  have  a  slight  anti-cyclonic  curve 
(see  Figure  79). 

2.  TEMPORAL  VARIATION 

To  better  illustrate  the  monthly  and  seasonal  variability  of  the  winds,  the  rate  of 
occurrence  of  the  six  patterns  was  computed  and  the  results  are  shown  in  Figure  80. 
Pattern  1  clearly  dominates  during  the  summer  when  the  weather  warms  up.  Its  maximum 
contribution  is  in  June  when  it  takes  up  almost  35%  of  the  data  set  and  it  disappears 
completely  in  January  and  February.  Pattern  2  is  more  likely  to  be  observed  during 
spring-like  and  fall-like  months.  Its  maximum  contribution  is  detected  in  May,  taking  up 
a  little  more  than  20%  of  the  data  set,  and  it  has  the  second  highest  percentage  among  the 
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Geostrophic  Wind  Pattern  :  6  (11.9658%) 
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six  patterns  in  this  particular  month.  This  pattern  declines  in  extreme  weather  conditions 
like  the  cold  weather  in  winter  and  the  warm  weather  in  the  summer.  Pattern  3  tends  to 
appear  mostly  during  the  second  six-month  period  of  the  year.  It  reaches  its  maximum  in 
August  and  becomes  the  pattern  most  likely  to  appear  during  that  time  of  the  year.  From 
July  to  mid-September  Pattern  3  is  the  most  dominant  pattern  among  the  six  patterns.  It 
tends  to  stay  at  15%-18%  during  the  first  six-month  period  of  the  year  and  never 
disappears  completely.  Pattern  4  comprises  the  highest  portion  of  the  data  set  during  the 
winter-like  and  summer-like  months.  It  is  the  most  dominant  pattern  in  February  with 
almost  35%  and  it  is  the  second  most  dominant  pattern  in  June  with  -28%.  Pattern  4  gets 
very  weak  during  the  spring  and  shows  an  oscillatory  pattern  of  contribution  during  the 
fall.  Pattern  5  contributes  significantly  during  the  spring  and  fall.  It  reaches  its  maximum 
percentages  in  fall-like  months  and  dominates  in  October  with  almost  40%  and  in  April 
with  -25%.  On  the  other  hand,  the  minimum  contribution  of  Pattern  5,  with  less  than 
10%,  is  observed  during  June.  Pattern  6  has  no  significant  seasonality  and  occurs 
occasionally  throughout  the  year.  It  contributes  the  least  among  the  six  patterns  from  May 
to  September  when  it  makes  up  only  about  5%  of  the  data  set.  With  a  little  less  than  30% 
contribution,  Pattern  6  is  the  pattern  most  likely  to  be  detected  in  March. 

The  time  series  (2000  to  2009)  of  the  six  wind  patterns  (shown  in  Figures  74-79) 
derived  by  SOM  is  shown  in  Figure  81,  and  the  evolution  of  the  duration  and  the 
frequency  of  the  same  patterns  is  shown  in  Figure  82  in  order  to  illustrate  the  inter-annual 
variability.  Furthermore,  Figure  83,  Figure  84  and  Figure  85  demonstrate  the 
contemporaneous  EAWR,  NAO,  and  ENSO  indices,  respectively,  along  with  the  times 
series  of  the  winds.  In  addition  to  these  figures,  Tables  10  and  11  provide  a  deeper 
understanding  of  the  relationship  between  the  wind  patterns  and  the  large-scale 
phenomena,  as  well  as  the  relationship  between  the  wind  patterns  and  the  surface  current 
patterns.  Figure  86  presents  the  anomalies  of  the  wind  components  for  each  month  (from 
2000  to  2009)  and  indicates  to  which  recurrent  patterns  of  the  surface  winds  (shown  in 
Figures  74-79)  these  anomalies  are  mapped  by  the  SOM  procedure. 

The  potential  influence  of  large-scale  teleconnection  indices  on  the  recurrent 
patterns  of  the  surface  winds  (shown  in  Figures  74-79)  are  evaluated  by  estimating  for 
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each  recurrent  pattern  the  value  of  the  climate  indices  (i.e.  EAWR,  NAO,  ENSO), 
acquired  for  the  different  months  mapped  towards  this  pattern.  This  process  enables 
determining  whether  a  recurrent  pattern  is  associated  to  either  positive  or  negative  value 
of  a  particular  index  significantly,  and  therefore  a  connection  between  the  surface  winds 
over  the  Black  Sea  and  large-scale  teleconnection  indices.  Figure  87  shows  for  each 
recurrent  pattern  of  the  surface  wind,  the  distribution  of  positive  and  negative  phases  of 
the  three  teleconnection  indices  during  the  months  mapped  towards  this  atmospheric 
pattern.  According  to  the  figure,  the  EAWR  index  and  the  ENSO  index  are  negative 
significantly  during  months  mapped  towards  pattern  2,  which  indicates  that  the 
occurrence  of  pattern  2  (The  Easterly  Winds  and  the  Anti-cyclonic  Curve  in  the  Western 
Shelf)  is  promoted  by  these  indices.  On  the  other  hand,  NAO  and  EAWR  indices  are 
significantly  positive  during  months  mapped  towards  pattern  4  (The  Northwesterly 
Winds  in  the  West  and  the  Westerly  Winds  in  the  East)  and  therefore  influence  the  wind 
anomaly. 
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Figure  80.  Monthly  percentage  of  the  six  wind  patterns 
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Geostrophic  Wind  Patterns  from  2000  to  2003 


Geostrophic  Wind  Patterns  from  2003  to  2006 


Geostrophic  Wind  Patterns  from  2006  to  2009 


Figure  8 1 .  The  inter-annual  variability  of  the  six  wind  patterns 


Year 


Figure  82.  Evolution  of  the  duration  and  the  frequency  of  the  six  patterns  of  the  winds 
shown  in  Figures  74-79.  For  instance,  the  SOM  Pattern  3  lasted  more  than  6  months  and 
dominated  approximately  2  months  between  years  2008  and  2009  with  more  than  50%. 
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Geostrophic  Wind  Patterns  from  2000  to  2009 


Figure  83.  Comparison  between  the  inter- annual  variability  of  the  six  wind  patterns  and 

the  EAWR  index 


Geostrophic  Wind  Patterns  from  2000  to  2009 


Figure  84.  Comparison  between  the  inter-annual  variability  of  the  six  wind  patterns  and 

the  NAO  index 
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Geostrophic  Wind  Patterns  from  2000  to  2009 


Figure  85.  Comparison  between  the  inter- annual  variability  of  the  six  wind  patterns  and 

the  ENSO  index 


Table  9.  The  percentages  of  the  six  wind  patterns  with  respect  to  the  large-scale 
phenomena.  For  example,  when  the  NAO  index  is  positive,  Pattern  3  has 
the  highest  possibility  to  form  among  the  six  patterns  at  32%.  On  the  other 
hand,  when  the  EAWR  index  is  positive,  Pattern  2  has  the  smallest 
possibility  of  occurring  at  7%. 


GEOS.  WIND 

(+)  NAO 

(-)  NAO 

(+)  EAWR 

(-)EAWR 

(+)ENSO 

(-)  ENSO 

Pattern  1 

11  % 

12% 

10% 

13  % 

11.5  % 

12.5  % 

Pattern  2 

11  % 

12% 

7  % 

15  % 

10% 

14% 

Pattern  3 

32  % 

18  % 

26  % 

23  % 

23  % 

26  % 

Pattern  4 

21  % 

13  % 

24% 

11  % 

19  % 

14% 

Pattern  5 

19  % 

25  % 

21  % 

24% 

25  % 

21  % 

Pattern  6 

6  % 

19% 

12% 

14% 

11.5  % 

12.5  % 
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Table  10.  The  percentages  of  the  six  wind  patterns  with  respect  to  the  six  patterns  of 
the  surface  currents.  For  example,  when  Pattern  3  of  the  winds  occurs, 
Pattern  2  of  the  surface  currents  has  the  highest  possibility  to  form  at 
5.54%.  On  the  other  hand,  when  Pattern  1  of  the  winds  occurs,  Pattern  4  has 
the  smallest  possibility  of  occurring  at  0.43%. 


GEOS.  WI> 

ID 

1 

2 

3 

4 

5 

6 

TOTAL 

GEOS.  VELOCTY 

1 

0.43  % 

3.2  % 

3.62  % 

3.62  % 

2.99  % 

3.2  % 

17.06  % 

2 

2.56  % 

1.92% 

5.54  % 

2.56  % 

4.9  % 

2.35  % 

19.83  % 

3 

3.62  % 

1.49  % 

3.2% 

1.28  % 

2.99  % 

1.07  % 

13.65  % 

4 

0.43  % 

2.35  % 

4.69  % 

3.2% 

4.69  % 

4.48  % 

19.83  % 

5 

2.77  % 

1.92% 

4.26  % 

2.56  % 

4.05  % 

1.49  % 

17.06  % 

6 

2.35  % 

0.85  % 

2.99  % 

3.2% 

2.56  % 

0.64  % 

12.57  % 

TOTAL 

12.15  % 

11.73  % 

24.31  % 

16.42  % 

22.17  % 

13.22  % 

100% 
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Figure  86.  Monthly  anomalies  of  the  components  of  zonal  and  meridional  winds.  The  six 
wind  patterns  provided  by  the  SOM  analysis  are  distinguished  by  the  large  circles.  The 
smaller  circles,  on  the  other  hand,  depict  the  basin  scale  averages  for  all  the  months  from 

January  2000  to  January  2009. 
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Figure  87.  Distribution  for  each  pattern  of  the  surface  currents  between  positive  (red)  and 
negative  (blue)  phases  of  the  teleconnection  indices.  The  areas  with  darker  colors 
represent  the  strong  phases  for  which  the  index’s  absolute  value  is  greater  than  a  standard 
deviation.  Monthly  average  values  of  the  climate  are  obtained  through  the  NOAA 
Climate  Prediction  Center  (Available  at  http://www.cpc.ncep.noaa.gov/data/teledoc/ 

telecontents .  shtml) . 
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VII.  CHLOROPHYLL 


Monthly  mean  climatological  chlorophyll-a  concentration  is  plotted 
logarithmically  and  these  results  are  displayed  in  Figures  88  to  99.  Additionally,  Figure 
100  presents  the  climatological  average  of  the  chlorophyll-a  concentration  from  1999  to 
2009.  In  January,  the  concentration  in  most  of  the  open  basin  stays  around  1  mg/m  . 
Higher  values  up  to  6  mg/m  are  observed  in  the  coastal  region  of  the  northwestern  shelf 
where  the  Danube,  Prut,  and  Dniester  rivers  flow  to  the  sea.  Likewise,  a  concentration 
level  up  to  5-6  mg/m  is  evident  in  the  north  region  of  the  Turkish  coast  where  the 
Kizilirmak  and  Yesilirmak  rivers  flow.  At  the  Crimean  Peninsula  coasts,  the 
concentration  slightly  decreases  and  values  stay  between  0.5  and  1  mg/m  (see  Figure 
88).  The  eastern  shelf  of  the  basin  maintains  a  steady  concentration  level  in  February,  but 
the  values  at  the  Western  shelf  generally  drop  slightly  except  in  the  coastal  region  at  the 
northwest  where  the  river  outflows  take  place.  In  this  particular  region,  higher  values  up 
to  5  mg/m  are  detected  and  the  area  of  this  higher  concentration  spreads  out  compared  to 
January.  Meanwhile,  high  concentration  levels  at  the  Turkish  coast  almost  completely 
disappear  (see  Figure  89).  Figure  90  shows  the  climatological  average  of  the  chlorophyll  - 
a  concentration  for  March.  The  chlorophyll-a  concentration  detected  is  approximately  0.5 
mg/m  basin-wide.  Small  regions  in  the  open  basin  reach  a  concentration  level  up  to  1 
mg/m  .  The  northwest  region,  which  is  influenced  by  the  river  outflows,  exhibits  high 
concentrations  up  to  4  mg/m  off  the  coast  and  up  to  10  mg/m  right  along  the  coast.  In 
April,  in  the  region  south  of  the  western  shelf  and  in  the  central  eastern  shelf  region, 
concentration  levels  of  1  mg/m  are  detected.  The  river  outflows  affect  a  larger  area  of  the 
northwestern  corner  of  the  basin,  but  the  highest  values  decrease  slightly  (they  are 
detected  at  6  mg/m3)  (see  Figure  91).  In  May,  the  concentration  levels  decrease 
dramatically  and  remain  at  approximately  0.2  mg/m  at  the  eastern  shelf  and  at  the  central 
western  shelf.  The  river  outflow  region  at  the  northwestern  comer  of  the  basin  elongates 
to  the  south  and  reaches  up  to  the  Bulgarian  coast  (see  Figure  92).  Figure  93  shows  the 
climatological  average  of  the  chlorophyll-a  concentration  during  June.  It  is  very  similar 
to  May  with  the  exception  of  the  large  area  at  the  northwestern  corner  that  is  influenced 
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by  the  river  outflows.  During  June,  higher  concentration  levels  (up  to  14  mg/m  )  are 
detected  in  that  region.  The  chlorophyll-a  concentration  levels  are  almost  identical  in 
July,  August  and  September.  The  concentration  level  at  the  northwestern  comer  begins  to 
shrink,  yet  the  maximum  chlorophyll-a  level  in  this  particular  region  is  still  very  high  as 
it  remains  around  10  mg/m  (see  Figure  94,  Figure  95,  and  Figure  96,  respectively). 
Figure  97  presents  the  climatological  average  of  the  chlorophyll-a  concentration  for 
October.  The  uniform  structure  of  the  basin  dissolves  and  a  chaotic  distribution  in  terms 
of  the  concentration  levels  takes  place.  Both  in  November  and  December,  the  chaotic 
distribution  of  the  chlorophyll-a  concentration  increases  in  intensity  and  scale.  The  river 
outflow  region  at  the  northwestern  corner  elongates  and  reaches  up  to  the  Bosphorus.  The 
central  coast  of  Turkey  exhibits  higher  than  average  values  as  well,  especially  where  the 
rivers  flow  to  the  Black  Sea  (see  Figure  98  and  Figure  99,  respectively).  Finally  the 
climatological  average  of  the  chlorophyll-a  concentration  from  1999  to  2009  is  shown  in 
Figure  100. 
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Figure  88.  The  climatological  average  of  the  chlorophyll-a  concentration  for  January 
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Figure  89.  The  climatological  average  of  the  chlorophyll-a  concentration  for  February 
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Figure  90.  The  climatological  average  of  the  chlorophyll-a  concentration  for  March 
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91.  The  climatological  average  of  the  chlorophyll-a  concentration  for  April 
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Figure  92.  The  climatological  average  of  the  chlorophyll-a  concentration  for  May 
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Figure  93.  The  climatological  average  of  the  chlorophyll-a  concentration  for  June 


Figure  94.  The  climatological  average  of  the  chlorophyll-a  concentration  for  July 
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Figure  96.  The  climatological  average  of  the  chlorophyll-a  concentration  for  September 
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Figure  95.  The  climatological  average  of  the  chlorophyll-a  concentration  for  August 
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97.  The  climatological  average  of  the  chlorophyll-a  concentration  for  Oct 
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Figure  98.  The  climatological  average  of  the  chlorophyll-a  concentration  for  November 
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Figure  99.  The  climatological  average  of  the  chlorophyll-a  concentration  for  December 
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Figure  100.  The  climatological  average  of  the  chlorophyll-a  concentration  from  2000  to 

2009 
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Figure  101  shows  the  horizontally  averaged  East  and  West  relative  vortices,  with 
the  chlorophyll-a  concentration  superimposed.  As  seen  in  the  figure,  Pattern  6  is 
characterized  by  low  chlorophyll-a  concentration  levels.  Aside  from  this  relation,  there  is 
no  indication  of  any  influence  of  the  vorticity  of  recurrent  surface  patterns  on 
chlorophyll-a  concentrations.  Likewise,  Figure  103  shows  monthly  anomalies  in  the 
components  of  zonal  and  meridional  winds,  with  the  chlorophyll-a  concentration 
superimposed.  Bahamon  et  al.  (2010)  state  that  there  is  no  significant  correlation  between 
the  wind  and  the  chlorophyll-a  level,  and  this  research  confirms  that  result.  There  is  no 
significant  separation  of  the  recurrent  patterns  of  the  wind  with  chlorophyll-a 
concentration  superimposed. 

Figure  102  and  Table  11  are  provided  to  allow  a  deeper  understanding  of  the 
correlation  between  the  recurrent  patterns  in  relation  to  both  spring  and  fall  blooms  of 
chlorophyll-a  concentrations.  According  to  the  figure  and  table,  Pattern  4  (Cyclonic  RIM 
Current  and  Cyclonic  Batumi  Eddy)  dominates  most  of  the  fall  blooms.  Apart  from  that, 
there  is  no  indication  of  any  correlation  concerning  the  spring  blooms.  Similarly,  Figure 
104  and  Table  12  show  the  same  kind  of  relation  between  the  recurrent  patterns  of  the 
surface  wind  and  the  chlorophyll-a  concentration.  But  this  time,  none  of  the  recurrent 
patterns  stand  out  and  the  distribution  is  random. 

Chu  et  al.  (2005)  used  the  optimal  spectral  decomposition  (OSD)  to  reconstruct 
seasonal  variability  of  the  Black  Sea’s  horizontally  averaged  chlorophyll-a  concentration 
from  data  collected  during  the  NATO  SfP-971818  Black  Sea  Project  in  1980-1995. 
According  to  their  study,  there  are  two  blooms  in  the  Black  Sea.  The  first  bloom  is  the 
winter/spring  bloom  which  takes  place  in  February  and  March,  and  terminates  at  the  end 
of  April.  The  second  bloom  is  the  fall  bloom  which  occurs  due  to  an  increase  in  the 
surface  chlorophyll-a  concentration  during  late  September.  The  fall  bloom  reaches  its 
maximum  in  October  and  has  half  the  intensity  of  the  winter/spring  bloom.  Chu  et  al.’s 
finding  contradicts  the  results  of  this  research,  which  determined  that  each  year  the  fall 
bloom  is  significantly  more  intense  than  the  winter/spring  bloom.  These  varying  findings 
can  be  explained  either  by  the  difference  in  the  data  sets  used  (in  situ  data  versus  satellite 
data)  or  the  difference  in  the  time  periods  during  which  each  study  occured.  Figure  105 

117 


shows  the  two-bloom  structure  of  the  surface  chlorophyll-a  concentration  Chu  et  al. 
discovered,  and  Figure  106  shows  the  distribution  of  the  observational  stations  they  used 
during  their  study. 
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Figure  101.  (a)  Horizontally  averaged  East  and  West  relative  vorticities,  with  (b) 
chlorophyll-a  concentration  superimposed.  Small  circles  represent  the  monthly  anomalies 
of  the  basin  scale  averaged  zonal  component  of  the  vorticity  of  the  surface  currents  from 
2000  to  2009.  Larger  circles  indicate  the  recurrent  spatial  patterns  identified  by  the  SOM 
analysis  that  are  illustrated  in  Figures  46-51.  The  basin  average  chlorophyll-a  anomaly 

is  superimposed. 
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Figure  102.  Temporal  evolution  of  the  chlorophyll-a  concentrations  from  1999  -  2009.  In 
order  to  assess  the  correlation  with  the  surface  currents,  the  bars  represent  the 
contribution  of  the  recurrent  patterns  for  each  three-month  period. 
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Table  11.  Distribution  of  the  occurrence  rates  of  the  surface  current  patterns  with 
respect  to  the  spring  and  fall  blooms  of  the  chlorophyll-a  concentrations 


SURFACE  CURR1 

ENT  PATTERNS 

YEAR 

BLOOM 

1 

2 

3 

4 

5 

6 

2000 

Spring 

20% 

- 

- 

- 

80% 

- 

Fall 

40% 

- 

15  % 

5  % 

40% 

- 

2001 

Spring 

5  % 

5  % 

65  % 

- 

- 

25  % 

Fall 

20% 

- 

- 

65  % 

15  % 

- 

2002 

Spring 

- 

100% 

- 

- 

- 

- 

Fall 

- 

- 

- 

100% 

- 

- 

2003 

Spring 

- 

70% 

5  % 

25  % 

- 

- 

Fall 

25  % 

15  % 

- 

60% 

- 

- 

2004 

Spring 

- 

65  % 

- 

35  % 

- 

- 

Fall 

20% 

- 

- 

- 

80% 

- 

2005 

Spring 

- 

- 

95  % 

- 

- 

5  % 

Fall 

- 

5  % 

- 

95  % 

- 

- 

2006 

Spring 

- 

55  % 

25  % 

- 

5  % 

15  % 

Fall 

100% 

- 

- 

- 

- 

- 

2007 

Spring 

40% 

50% 

5  % 

- 

5  % 

- 

Fall 

30% 

- 

- 

70% 

- 

- 

2008 

Spring 

- 

- 

45  % 

- 

- 

55  % 

Fall 
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Figure  103.  In  the  right  panel,  small  circles  represent  the  monthly  anomalies  of  the  basin 
scale  averaged  meridional  and  zonal  components  of  the  winds  from  2000  to  2009.  Larger 
circles  indicate  the  recurrent  spatial  patterns  identified  by  the  SOM  analysis  that  are 
illustrated  in  Figures  74  -  79.  The  basin  average  chlorophyll  anomaly  is  superimposed. 
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Year 


Figure  104.  Temporal  evolution  of  the  chlorophyll-a  concentrations  from  1999  -  2009.  In 
order  to  assess  the  correlation  with  the  surface  winds,  the  bars  represent  the  contribution 
of  the  recurrent  patterns  for  each  three-month  period. 


Table  12.  Distribution  of  the  occurrence  rates  of  the  recurrent  patterns  with  respect  to 
the  spring  and  fall  blooms  of  the  chlorophyll-a  concentrations. 


RE  Cl 

LJRRENT  PATTERNS 

YEAR 

BLOOM 

1 

2 

3 

4 

5 

6 

2000 

Spring 

20% 

5  % 

35  % 

20% 

10% 

10% 

Fall 

- 

15  % 

15  % 

5  % 

25  % 

40% 

2001 

Spring 

15  % 

5  % 

15  % 

15  % 

45  % 

5  % 

Fall 

- 

5  % 

35  % 

25  % 

15  % 

20% 

2002 

Spring 

10% 

10% 

30% 

20% 

30% 

- 

Fall 

- 

- 

50% 

5  % 

25  % 

20% 

2003 

Spring 

40% 

15  % 

25  % 

- 

15  % 

5  % 

Fall 

- 

15  % 

5  % 

40% 

15  % 

25  % 

2004 

Spring 

5  % 

5% 

45  % 

5  % 

20% 

20% 

Fall 

- 

15  % 

- 

55  % 

30% 

- 

2005 

Spring 

45  % 

- 

5  % 

20% 

15  % 

15  % 

Fall 

- 

5  % 

35  % 

10% 

20% 

30% 

2006 

Spring 

25  % 

5  % 

40% 

25  % 

5  % 

- 

Fall 

5  % 

5  % 

15  % 

25  % 

15  % 

35  % 

2007 

Spring 

20% 

45  % 

- 

15  % 

20% 

- 

Fall 

- 

15  % 

35  % 

15  % 

30% 

5  % 

2008 

Spring 

40% 

20% 

5  % 

25  % 

5  % 

5  % 

Fall 

5  % 

25  % 

5  % 

5  % 

45  % 

20% 
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Figure  105.  Sensitivity  of  the  temporal  interruption  on  chlorophyll-a  concentration  using 
(1)  all  data,  (2)  all  data  except  April,  and  (3)  all  data  except  August  (from  Chu  et  al. 

2005). 
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Figure  106.  Distribution  ofobservational  stations:  (a)  winter,  (b)  spring,  (c)  summer, and 

(d)  fall  (from  Chu  et  al.  2005). 
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VIII.  CONCLUSIONS 


Surface  geostrophic  velocity  data  from  1999  to  2009  were  obtained  from  AVISO 
to  provide  data  on  the  surface  currents.  Their  monthly  means  show  the  seasonal  variation 
of  the  surface  currents.  Six  patterns  were  found  after  conducting  the  SOM.  The  inter¬ 
annual  variability  of  the  surface  currents  was  also  obtained.  The  first  pattern  (Sevastopol 
Cyclonic  and  Batumi  Dipole  Eddies)  represents  more  than  20%  of  the  whole  data  set.  In 
this  pattern,  the  Batumi  dipole  eddies  are  located  at  the  southeastern  corner  of  the  basin. 
At  the  northwestern  corner  of  the  basin,  the  cyclonic  Sevastopol  eddy  is  evident.  The 
general  circulation  is  comprised  of  two  main  gyres.  The  western  main  gyre  is  anti- 
cyclonic  and  weak.  The  eastern  main  gyre,  on  the  other  hand,  is  cyclonic  and  strong.  The 
northerly  boundary  current  forms  in  the  west.  No  eddy  formation  is  observed  in  the  open 
parts  of  the  basin.  More  than  16%  of  the  data  set  is  represented  by  the  second  pattern 
(Cyclonic  RIM  Current  and  Anti-cyclonic  Batumi  Eddy).  The  main  feature  of  this  pattern 
is  the  strong,  anti-cyclonic  Batumi  eddy.  Apart  from  the  Batumi  eddy,  no  other  eddy 
structure,  including  the  major  Sevastopol  eddy,  is  observed.  The  general  circulation  is 
formed  by  the  strong  and  cyclonic  RIM  current.  The  open  parts  of  the  basin  are  relatively 
calm  as  opposed  to  the  coastal  regions.  Therefore,  this  particular  pattern  is  mostly 
dominated  by  the  RIM  current.  The  third  pattern  (Anti-cyclonic  Sevastopol  and  Batumi 
Eddies)  also  represents  almost  16%  of  the  whole  data  set.  The  major  eddies  of  the  basin 
are  the  Batumi  and  Sevastopol  eddies.  Both  of  these  major  eddies  are  strong  and  anti- 
cyclonic.  The  RIM  current  does  not  appear  which  means  that  the  general  circulation  is 
formed  by  the  western  and  eastern  main  gyres.  In  the  northeastern  corner,  the  weak 
Caucasus  eddy  exists,  while  the  strong  northerly  boundary  current  is  present  in  the  west. 
The  fourth  pattern  (Cyclonic  RIM  Current  and  Cyclonic  Batumi  Eddy)  comprises  more 
than  20%  of  the  whole  representation  of  the  data  set.  In  this  pattern  the  Black  Sea  basin  is 
almost  entirely  dominated  by  the  strong  cyclonic  RIM  current.  The  open  parts  of  the  sea 
stay  relatively  calm  and  the  very  weak  cyclonic  Batumi  eddy  is  detected  but  it  is  nearly 
absorbed  by  the  RIM  current.  The  fifth  pattern  (Anti-cyclonic  RIM  Current  and  Batumi 
Dipole  Eddies)  represents  15%  of  the  whole  data  set  and  it  demonstrates  a  basin-wide, 
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chaotic  environment.  At  the  southeastern  corner  the  Batumi  dipole  eddies  occur.  The 
weak,  anti-cyclonic  Crimea  eddy  is  observed  in  the  north,  whereas  the  Sevastopol  eddy 
disappears.  The  general  circulation  is  formed  by  the  strong,  anti-cyclonic  RIM  current. 
More  than  10%  of  the  total  representation  for  the  data  set  is  covered  by  the  sixth  pattern 
(Anti-cyclonic  RIM  Current  and  Multi  Eddies).  The  southeastern  corner  of  the  basin  is 
entirely  dominated  by  the  strong,  anti-cyclonic  Batumi  eddy.  At  the  north  of  the  Batumi 
eddy,  the  cyclonic  Suchumi  eddy  is  detected.  The  weak,  anti-cyclonic  Kerch  eddy  is 
observed  in  the  north.  At  the  northwestern  corner  the  strong  anti-cyclonic  Sevastopol 
eddy  forms.  The  general  circulation  is  comprised  of  the  strong,  anti-cyclonic  RIM 
current. 

The  rate  of  occurrence  of  the  six  patterns  was  computed.  Pattern  1  is  one  of  the 
two  patterns  that  dominate  during  winter-like  and  fall-like  months.  Its  maximum 
contribution  is  in  January  with  -45%  and  it  disappears  completely  in  lune.  Pattern  2  is 
more  likely  to  be  observed  during  the  spring.  Its  maximum  contribution  is  detected  in 
April  with  more  than  50%.  It  declines  in  winter-like  months,  and  it  completely  disappears 
in  June.  Pattern  3  tends  to  appear  towards  the  end  of  the  spring  and  the  beginning  of  the 
summer.  It  reaches  its  maximum  in  May  and  becomes  the  pattern  most  likely  to  appear 
during  that  time  of  the  year.  The  first  three  months  of  the  year  are  dominated  by  the 
Pattern  4.  It  tends  to  appear  strongly  during  cold  periods,  whereas  during  warm  periods  it 
exhibits  a  weak  contribution.  Pattern  4  disappears  in  May  and  doesn’t  appear  again  until 
September.  Despite  the  fact  that  Pattern  5  shows  maximum  percentages  in  fall-like 
months  and  dominates  in  October,  it  has  no  significant  seasonality  and  occurs 
occasionally  throughout  the  year.  Pattern  6  completely  disappears  from  November  to 
April,  and  appears  only  five  months  during  the  year.  It  dominates  the  period  from  June  to 
September  and  reaches  its  maximum  in  July.  With  more  than  40%  contribution,  it 
becomes  the  pattern  most  likely  to  be  detected  in  July.  The  EAWR  index  and  the  ENSO 
index  are  detected  significantly  negative  during  months  mapped  towards  pattern  5,  which 
indicates  that  the  occurrence  of  pattern  5  is  promoted  by  these  indices.  On  the  other  hand, 
ENSO  index  is  detected  significantly  positive  during  months  mapped  towards  pattern  2 
and  therefore  influence  the  surface  current  structure. 
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The  wind  data  from  1992  to  2012  were  obtained  from  QuikSCAT.  This  data  is 
representative  of  the  winds  near  the  surface  of  the  Black  Sea.  The  data  provides  monthly 
means  that  show  the  seasonal  variation  of  the  surface  currents.  After  conducting  the 
SOM,  six  patterns  were  found.  Five  out  of  these  six  recurrent  patterns  perfectly  match 
with  the  results  of  the  study  of  Arthur  et  al.  (2012).  Additionally,  the  inter- annual 
variability  of  the  winds  was  obtained.  This  research  focused  on  the  time  frame  from  1999 
to  2009. 

The  first  pattern  (The  Northerly  Winds  and  the  Cyclonic  Curve  in  the  Western 
Shelf)  represents  more  than  13%  of  the  whole  data  set.  In  this  pattern,  westerly  winds  are 
observed  north  of  the  western  shelf.  The  winds  curve  anti-cyclonically  and  become 
northeasterly  south  of  the  western  shelf.  At  the  eastern  self,  on  the  other  hand,  the  winds 
blow  from  the  northwest  with  a  slight  cyclonic  curve  south  of  the  eastern  shelf.  More 
than  13%  of  the  total  representation  for  the  data  set  is  covered  by  the  second  pattern  (The 
Easterly  Winds  and  the  Anti-cyconic  Curve  in  the  Western  Shelf).  The  main  feature  of 
this  pattern  is  the  anti-cyclonic  curve  of  the  winds  in  the  western  shelf.  The  winds  blow 
from  the  east,  starting  at  the  eastern  shelf,  and  from  the  central  Black  Sea  region.  The 
winds  are  stronger  to  the  west  and  they  are  strongest  north  of  the  western  shelf.  The  third 
pattern  (The  Dominance  of  the  Strong  Northerly  Winds)  represents  23%  of  the  whole 
data  set.  The  strong  northerly  winds  almost  completely  dominate  the  Black  Sea  region. 
The  winds  lose  their  strength  to  the  south  of  the  eastern  shelf  and  they  exhibit  a  cyclonic 
curve  in  this  region.  The  fourth  pattern  (The  Northerly  Winds  in  the  West  and  the 
Westerly  Winds  in  the  East)  comprises  more  than  16%  of  the  whole  representation  of  the 
data  set.  The  western  and  eastern  shelves  exhibit  two  different  wind  characteristics.  The 
strong  northwesterly  winds  occur  in  the  western  shelf  with  a  little  anti-cyclonic  curve  to 
the  south,  whereas,  starting  from  the  central  Black  Sea,  the  winds  blow  from  the  west  and 
the  whole  eastern  shelf  is  dominated  by  the  relatively  weak  easterly  winds.  The  fifth 
pattern  (The  Dominance  of  the  Strong  Northeasterly  Winds)  represents  more  than  21%  of 
the  whole  data  set.  The  northeasterly  winds  dominate  the  whole  Black  Sea  region  with 
the  exception  of  the  southeastern  corner  of  the  basin.  The  winds  become  weaker  to  the 
southeast  and  almost  disappear  completely  close  to  the  coast.  Almost  12%  of  the  total 
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representation  for  the  data  set  is  covered  by  the  sixth  pattern  (The  Strong  Southwesterly 
Winds).  The  southwesterly  winds  dominate  the  whole  Black  Sea  region  with  the 
exception  of  the  southeastern  comer  of  the  basin.  The  winds  blowing  from  the  southeast 
become  weaker  and  they  have  a  slight  anti-cyclonic  curve. 

The  rate  of  occurrence  of  the  six  patterns  is  computed.  Pattern  1  clearly  dominates 
during  the  summer  when  the  weather  warms  up.  Its  maximum  contribution  is  in  June  and 
it  disappears  completely  in  January  and  February.  Pattern  2  is  more  likely  to  be  observed 
during  spring-like  and  fall-like  months.  Its  maximum  contribution  is  detected  in  May. 
This  pattern  declines  in  extreme  weather  conditions  like  the  cold  weather  in  winter  and 
the  warm  weather  in  the  summer.  Pattern  3  tends  to  appear  mostly  during  the  second  six- 
month  period  of  the  year.  It  reaches  its  maximum  in  August  and  becomes  the  pattern 
most  likely  to  appear  during  that  time  of  the  year.  It  never  disappears  completely.  Pattern 
4  comprises  the  highest  portion  of  the  data  set  during  the  winter-like  and  summer-like 
months.  It  is  the  most  dominant  pattern  in  February.  Pattern  4  gets  very  weak  during  the 
spring  and  shows  an  oscillatory  pattern  of  contribution  during  the  fall.  Pattern  5 
contributes  significantly  during  the  spring  and  fall.  It  reaches  its  maximum  percentages  in 
fall-like  months  and  dominates  in  October.  Pattern  6  has  no  significant  seasonality  and 
occurs  occasionally  throughout  the  year.  It  contributes  the  least  among  the  six  patterns 
from  May  to  September.  Pattern  6  is  the  pattern  most  likely  to  be  detected  in  March.  The 
EAWR  index  and  the  ENSO  index  are  detected  significantly  negative  during  months 
mapped  towards  pattern  2,  which  indicates  that  the  occurrence  of  pattern  is  promoted  by 
these  indices.  On  the  other  hand,  NAO  and  EAWR  indices  are  detected  significantly 
positive  during  months  mapped  towards  pattern  4  and  therefore  influence  the  wind 
anomaly. 

According  to  the  horizontally  averaged  East  and  West  relative  vortices,  with  the 
chlorophyll-a  concentration  superimposed,  Pattern  6  is  characterized  by  low  chlorophyll- 
a  concentration  levels.  Aside  from  this  relation,  there  is  no  indication  of  any  influence  of 
the  vorticity  of  recurrent  surface  patterns  on  chlorophyll-a  concentration.  On  the  other 
hand,  the  monthly  anomalies  in  the  components  of  zonal  and  meridional  winds,  with  the 
chlorophyll-a  concentration  superimposed  show  no  significant  separation,  hence  no 
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significant  correlation  between  the  wind  and  the  chlorophyll-a  level.  This  finding 
confirms  the  study  of  Bahamon  et  al.  (2010).  Pattern  4  of  the  surface  current  data  set 
dominates  most  of  the  fall  blooms.  Apart  from  that,  there  is  no  indication  of  any 
correlation  concerning  the  spring  blooms.  On  the  other  hand,  none  of  the  recurrent 
patterns  of  the  wind  data  set  stands  out;  therefore  the  distribution  is  determined  as 
random. 

Chu  et  al.  (2005)  found  two  blooms  in  the  Black  Sea  using  in-situ  measurements 
in  1980-1995.  The  first  bloom  is  the  winter/spring  bloom  which  takes  place  in  February 
and  March,  and  terminates  at  the  end  of  April.  The  second  bloom  is  the  fall  bloom  which 
occurs  due  to  an  increase  in  the  surface  chlorophyll-a  concentration  during  late 
September.  The  fall  bloom  reaches  its  maximum  in  October  and  has  half  the  intensity  of 
the  winter/spring  bloom.  Our  study  shows  the  change  of  the  bi-modal  characteristics  in 
1999-2009  with  the  fall  bloom  being  more  significant  than  the  spring  bloom.  The  surface 
circulation  Pattern  4  (cyclonic  RIM  current  and  Batumi  eddy)  is  associated  with  the 
occurrence  of  the  fall  bloom.  Evident  connection  of  negative  NAO  and  negative  ENSO  to 
the  Pattern  4  circulation  implies  the  large-scale  atmospheric  effect. 
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